Robotics  Research 
Ibchnical  Report 


Partitioning  Arrangements  of  Lines:  I.  An 
Efficient  Deterministic  Algorithm 

by 

Pankaj  K.  Agarwal 

Technical  Report  No.  446 

Robotics  Report  No.  197 

May,  1989 


♦• 


1 1 


\, 


o 


CO 

c 

e 

0) 

en 

c 

\<  (0 
u 
•n  u 
(0  (0 
J^  • 

C     CPM 

(C  c 

Plj  -H    «• 

C   CQ 

-  O   CU 

-H  -H    C 

S  -H  iH 
U  4J 
(0    i-l  M-l 
Cn  (0    O 


o 


New  York  University 
int  Institute  of  Mathematical  Sciences 

Computer  Science  Division 

25 1  Mercer  Street  New  York,  N.Y  1 00 1 2 


:>: 


it.  * --^p 


Partitioning  Arrangements  of  Lines:   I.  An 
Efficient  Deterministic  Algorithm 

by 

Panlcaj  K.  Agarwal 


Technical  Report  No.  446 

Robotics  Report  No.  197 

May,  1989 


New  York  University 

Dept.  of  Computer  Science 

Courant  Institute  of  Mathematical  Sciences 

251  Mercer  Street 

New  York,  New  York   10012 


Work  on  this  paper  has  been  supported  by  Office  of  Naval  Research  Grant  N00014-87-K- 
0129  National  Science  Foundation  CER  Grant  DCR-83-20085,  National  Science  Foundntinn 
Grant  subcontract  CMU-406349-55586,  and  by  grants  from  the  Digital  Equipment  Corpcr  >- 
tion  and  the  IBM  Corporation. 


Partitioning  Arrangements  of  Lines:  I.  An 
Efficient  Deterministic  Algorithm 

Pankaj    K.  AgaswaP 

Courant  Institute  of  Mathematical  Sciences 

New  York  University 

NY  10012 


ABSTRACT 

In  this  paper  we  consider  the  following  problem:  Given  a  set  C  oi  n  lines  in 
the  plane,  partition  the  plane  into  O(r^)  triangles  so  that  no  triangle  meets  more 
than  0{j)  lines  of  £.  We  present  a  deterministic  algorithm  for  this  problem  with 
0(rarlognlog'^  r)  running  time,  where  u;  is  a  constant  <  3.33. 
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Introduction 


1      Introduction 

In  the  last  few  years  several  randomized  divide  and  conquer  algorithms  for  a  variety  of 
geometric  problems  have  been  developed  {[AS],  [Cla],  [Clb],  [Clc],  [CS],  [CTV],  [EGS], 
[GOS]  etc.),  using  e-nets  [HW]  or  the  random  sampUng  technique  of  [Clb]  (see  also  [RSa], 
[RSb]).  The  e-net  theory  shows  that  for  a  given  set  X  of  n  objects  and  a  set  7^  C  2^  of 
ranges  with  finite  Vapmk-Chervonenkis  dimension  (see  [HW]  for  definition),  there  exists 
a  subset  A''  C  X  of  size  r,  for  any  r  >  0,  such  that  if  r  G  7^  and  t  f]  N  =  %,  then  \t\ 
=  O(-logr).  Moreover,  the  theory  asserts  that  a  random  saxnple  A''  of  size  r  will  be  an 
e-net  with  high  probability.  A  similar  result  has  been  independently  obtained  by  Clarkson 
[Clb],  [Clc].  It  states  roughly  that  if  we  draw  a  random  sample  A''  of  size  r  from  our 
set  X  and  partition  the  underlying  geometric  space  into  cells  with  the  property  that  (i) 
each  cell  can  be  defined  in  terms  of  only  a  "constant"  number  of  elements  of  A'',  and  (ii) 
each  cell  does  not  meet  any  element  of  A'',  then,  with  high  probability,  no  cell  meets  more 
than  0(^  logr)  elements  of  X.  Random  sampling  techniques  allow  us  to  split  the  original 
problem  into  subproblems  of  small  size,  which  can  then  be  solved  recursively.  For  example, 
if  the  objects  are  "lines  in  the  plane"  and  em  element  of  7^  is  the  "set  of  lines  intersecting  a 
given  triangle",  then  these  results  imply  that,  for  a  given  set  £  of  n  lines  in  the  plane,  there 
exists  a  subset  N  C  C  oi  size  r  such  that  any  triangle  which  misses  all  lines  in  A''  intersects 
at  most  O(Mogr)  lines  of  £.  A  problem  involving  the  ^et  £  of  lines  can  then  be  spht  into 
subproblems  by  computing  the  arrangement  of  the  lines  in  N  and  by  triangulating  each 
face  of  the  axrangement;  each  resulting  triangle  A  induces  a  subproblem  involving  only 
O(Mogr)  lines  of  £  intersecting  A.  As  an  example  of  a  problem  that  benefits  from  this 
set-up,  consider  Hopcroft'a  problem  (see  [CSY]):  Given  m  points  and  n  Unes  in  the  plane, 
does  any  point  lie  on  any  Une?  Applying  the  above  divide  and  conquer  strategy  to  the  given 
lines,  and  partitioning  the  points  among  the  triangles,  each  subproblem  involves  only  the 
points  that  fall  within  some  triangle  and  the  line  cutting  it.  More  details  concerning  this 
and  other  problems  of  this  kind  are  given  in  the  second  paxt  of  this  paper  [Aga].  The  same 
approach  works  for  other  geometric  objects  in  the  plane  as  well  as  in  higher  dimensions. 

Unfortunately  the  random  sampling  technique  or  the  e-net  theory  only  proves  the 
existence  of  such  a  subset  N,  but  does  not  show  how  to  construct  it  deterministicaily  (and 
efficiently).  This  has  forced  ail  algorithms  based  on  these  techniques  to  use  randomization 
to  obtain  N  (the  only  randomized  step  in  most  of  these  algorithms  is  that  of  choosing 
a  good  sample  N).  The  randomization  is  justified  because  the  above  theories  show  that 
most  subsets  N  of  size  r  satisfy  the  desired  properties,  so  a  random  choice  of  such  a 
subset  will  succeed  with  high  probability.  Still,  a  major  open  problem  in  this  area  is  to 
find  an  efiicient  deterministic  algorithm  to  obtain  a  "good"  sample.  Recently  Chazelle 
and  Friedman  [CF]  gave  a  general  framework  which  unifies  the  results  of  Haussler-Welzl 
[HW]  and  Clarkson  [Clc],  and  yields  a  deterministic  algorithm  to  construct  a  good  sample. 
Although  their  algorithm  runs  in  polynomial  time,  it  is  not  very  efficient  (in  particular  its 
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time  complexity  is  too  high  an  overhead  for  most  of  the  apphcations).  This  has  motivated 
researchers  to  look  for  special  ca^es,  where  faster  deterministic  constructions  are  possible. 
Woeginger  [Wo]  gave  an  0{n)  algorithm  to  construct  an  e-net  for  a  very  special  case, 
where  the  objects  are  points  in  the  plane  and  the  ranges  are  half-planes.  A  much  more 
significant  result  in  this  direction  has  been  obtained  by  Matousek  [Maa],  who  showed  that, 
given  a  set  of  n  hnes  in  the  plane,  the  plane  can  be  partitioned  into  0(r'^)  triangles  in  time 
0{nr^  log^  r)  so  that  no  triangle  intersects  more  than  0{^)  Unes  (however,  these  triangles 
in  general  are  not  formed  by  an  arrangement  of  a  sample  subset  of  the  lines).  Note  that 
an  important  property  of  the  algorithm  is  that  each  triangle  it  produces  is  cut  by  only 
O(-)  lines,  rather  than  O(jlogr),  as  promised  by  the  probabihstic  techniques.  If  r  is 
constant,  then  this  algorithm  is  optimal,  and  can  be  used  to  remove  randomization  from 
most  of  the  "random- sampling"  based  algorithms  involving  lines  or  segments  ([AS],  [EGS], 
[GOS],  [HW]  etc.)  However,  it  is  not  efficient  if  r  is  large,  in  particular  if  r  =  n",  for 
some  0  <  a  <  1.  Let  us  be  more  specific  concerning  this  efficiency  issue.  Even  though 
Matousek's  algorithm  is  much  faster  than  that  of  [CF],  it  is  still  inefficient  when  compared 
to  the  following  "lower  bound".  Since  the  algorithm  produces  Oir"^)  triangles  and  each  is 
met  by  O(-)  Unes,  the  total  input  size  of  all  the  corresponding  subproblems  (that  is,  the 
total  number  of  line-triangle  crossings)  is  0{nr)  (one  can  show  that  the  bound  is  tight  in 
the  worst  case).  Thus  it  is  natural  to  seek  an  improved  algorithm  whose  complexity  is 
close  to  0{nr).  As  we  will  see  later,  this  improvement  does  malce  a  difference,  because  in 
many  applications  we  do  want  to  choose  r  to  be  n",  for  some  0  <  a  <  1. 

In  this  paper  we  achieve  such  an  improvement,  and  obtain  an  C)(nr  log  n  log'^  r)  deter- 
ministic algorithm  which,  given  a  set  £  of  n  Unes,  partitions  the  plane  into  0{r^)  trismgles, 
none  of  which  meets  more  than  0{j)  lines  (here  u  is  some  constant  <  3.33).  Then  we 
apply  this  algorithm  to  remove  randomization  and  to  improve  the  time  complexity  of  so- 
lutions to  several  problems.  A  common  characteristic  of  the  time  complexity  of  previously 
known  algorithms  for  these  problems  is  that  they  axe  worse  by  a  factor  of  0{n^),  for  any 
6  >  0  (with  a  constant  of  proportionality  depending  on  6),  than  the  worst  case  output 
size,  or  than  some  other  natural  meastire  of  complexity  that  one  can  associate  with  the 
problem.  For  example,  the  running  time  of  the  best  known  (randomized)  algorithm  to 
compute  incidences  between  m  points  and  n  Unes  is  0(m^/^~*n^/^''"^*  +  {m  +  n)  log  n),  for 
any  6  >  0  [EGS],  which  is  worse  by  a  factor  of  0((^)*)  than  the  maximimi  possible  num- 
ber of  such  incidences.  (A  sUghtly  faster  but  still  randomized  algorithm  for  this  problem 
has  been  given  in  [EGH*]  which  nms  in  0{m'^^^n'^^^  \og^  n  +  {m  +  n^/^)log^n)  expected 
time.  The  difference  between  these  two  algorithms  is  that  the  first  adgorithm  can  be  made 
deterministic  without  increasing  its  ninning  time,  using  Matousek's  algorithm,  while  the 
second  algorithm  is  not  yet  known  to  admit  a  similarly  efficient  determinization.)  Factors 
Uke  0{n^)  appear  in  the  bound  of  these  algorithms  due  to  the  fact  that  they  use  small 
values  of  r  (in  most  cases  constant  values).  Although  Matousek's  algorithm  makes  most  of 
these  algorithms  deterministic,  it  also  cannot  choose  a  large  value  of  r,  because  its  running 
time  is  quadratic  in  terms  of  r,  in  which  case  the  overhead  of  constructing  the  partition 
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dominates  substantially  the  time  complexity,  resulting  in  an  inefficient  algorithm.  On  the 
other  hand,  by  applying  our  algorithm,  we  can  use  a  sufficiently  large  value  of  r,  which 
allows  us  to  remove  the  0(n*)  factor  from  the  bounds.  Another  disadvantage  of  using  a 
small  value  of  r  is  that  the  algorithm  then  becomes  recursive  and  more  complex.  In  con- 
trast, by  choosing  an  appropriate  large  value  of  r  we  partition  the  plane  just  once,  because 
then  the  subproblems  axe  sufficiently  small,  ajid  can  be  solved  directly  by  other  means; 
this  malces  the  algorithms  much  simpler. 

In  a  companion  paper  [Aga]  we  describe  several  applications  of  our  partitioning  algo- 
rithm including  the  following  ones: 

(i)  Computing  incidences  between  lines  and  points:  Given  a  set  of  n  lines  and  a  set  of  m 
points  in  the  plane,  compute  how  many  lines  pass  through  each  point.  (Alternatively, 
compute  the  lines  passing  through  each  point,  or  just  determine  whether  any  line 
passes  through  any  point.)  Edelsbrunner  et  al.  [EGS]  have  given  a  randomized  algo- 
rithm for  this  problem  with  expected  running  time  0{m^^^~^n^^'^'*'^^  +  {m  +  n)  log  n), 
for  any  ^  >  0.  A  slightly  improved,  but  still  randomized,  algorithm  has  been  given  in 
[EGH*].  We  present  a  deterministic  algorithm  with  Oim^^^n^^^log^^^ nlog'^^^  ^  + 
{m  +  n)logn)  time  complexity.  Since  the  meiximum  number  of  incidences  between 
m  points  and  n  lines  is  Q{m?l^n^l^  -\-Tn-\-  n),  our  algorithm  is  close  to  optimal  in  the 
worst  case. 

(ii)  Computing  many  faces  in  an  arrangement  of  lines:  Given  a  set  of  n  lines  and  a  set  of 
m  points  in  the  plane,  compute  the  faces  in  the  arrangement  of  the  lines  containing 
the  given  points.  Edelsbrunner  et  al.  [EGS]  have  given  a  randomized  algorithm  for 
this  problem  with  expected  ninning  time  C>(m^/^~'^n^/^+^*  -|-  nlognlogm),  for  any 
S  >  0.  As  in  the  case  of  the  incidence  problem,  a  slightly  better  randomized  algorithm 
has  been  given  in  [EGH*].  We  present  a  deterministic  0{m^^^n'^^^  log^/^  n  log*^^^  -^  -|- 
nlogn)  edgorithm,  again  coming  close  to  optimal  in  the  worst  case  (see  [CEGSW]). 

(iii)  Computing  many  faces  in  an  arrangement  of  segments:  This  is  the  same  problem  as 
the  previous  one  except  that  now  we  have  a  collection  of  segments  instead  of  lines. 
The  previous  best  solution  is  by  Edelsbrunner  et  al.  [EGS],  which  is  randomized  and 
has  expected  running  time  0(m^/^~*n^/^+^*-hna(n)  log^  n  log  m),  for  any  6  >  0,  where 
a{n)  is  a  fimctional  inverse  of  Ackermann's  function.  We  present  a  deterministic 
algorithm  with  improved  time  complexity  0(m^/^n^/^  log  nlog'^^^''"^  -y-  +  nlog^  n). 

(iv)  Counting  segment  intersections:  We  give  a  deterministic  O(n''/^log*'^'''^'/^n)  algo- 
rithm to  coimt  the  number  of  intersections  in  a  given  collection  of  n  segments;  this 
is  an  improvement  over  Guibas  et  al.'s  algorithm  [GOS],  which  counts  the  number 
of  intersections  in  0{n*^^'^^)  randomized  expected  time,  for  any  S  >  0. 

(v)   Counting  and  reporting  red-blue  intersections:  Given  a  set  F,  of  n^  "red"  segments 
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and  another  set  Fj  of  rif,  "blue"  segments  in  the  plane,  count  the  number  of  intersec- 
tions between  Fr  and  F^,  or  report  all  of  them.  (In  this  problem,  we  need  to  ignore 
the  potentially  large  number  of  intersections  within  Fr  or  within  Ff,.)  The  previous 
best  solution  is  by  Agaxwal  and  Sharir  [AS],  which  reports  all  K  red-blue  intersec- 
tions deterministically  in  0((nrv^  +  n(,y^-f  A")  log  n)  time,  where  n  =  Ur  +  rik.  We 
give  a  deterministic  O(n''/^log^'^'''^^^^n)  algorithm  to  count  all  red-blue  intersections. 
It  can  report  all  K  red-blue  intersections  in  0{n^l^  log^'^"^^^/^  n  -|-  K)  time. 

(vi)  Implicit  point  location  problem:  Given  a  collection  of  m  points  and  a  collection  of 
(possibly  intersecting)  n  triangles  in  the  plane,  find  which  points  lie  in  the  union  of 
the  triangles.  This  turns  out  to  be  a  special  ca^e  of  a  general  problem  of  implicit  point 
location  of  planar  maps  formed  by  overlapping  figures.  We  present  a  deterministic 
algorithm  with  0(m^/^n^/^  log^^^  n  log'^''^  -^  +  {m  +  n)logn)  time  complexity. 

(vii)  Approximate  hadf-plane  range  query:  Given  a  set  5  of  n  points  in  the  plane  and  a 
parameter  (not  necessarily  constant)  e  >  0,  preprocess  them  so  that,  for  any  query 
line  £,  we  can  approximately  count  the  number  of  points  lying  above  i  with  an  error 
of  at  most  ±en.  We  give  an  algorithm  that  preprocesses  5,  in  time  0(-  lognlog'^  -), 
into  a  data  structure  of  size  0(  j),  so  that  a  query  can  be  answered  in  0(log  n)  time. 

(viii)  Constructing  spanning  trees  with  low  stabbing  number.  Given  a  set  5  of  n  points 
in  the  plane,  we  present  an  0(n''/^  log"'"''^  n)  algorithm  to  construct  a  family  of 
k  =  O(logn)  spanning  trees  7i,  . . .  ,7jfe  of  5  with  the  property  that,  for  any  Une 
£  there  is  tree  Ti,  such  that  £  intersects  at  most  0(\/n)  edges  of  7^.  Moreover,  with 
additional  0(n  log  n)  preprocessing  and  0{n)  space,  the  tree  7^  corresponding  to  a 
query  hne  £  can  be  determined  in  O(log  n)  time.  The  previous  best  known  cdgorithm 
is  by  Matousek  [Mab],  which  runs  in  0(n^/'' log^  n)  time,  and,  moreover,  produces  a 
stabbing  number  0{y/n  \og^  n)  instead  of  0{y/n). 

(ix)  Space  query  time  tradeoff  in  simplex  range  searching:  Given  a  set  5  of  n  points 
in  the  plane  and  a  parameter  n^'^^  <  rn  <  n^~'',  for  some  constants  eo,ei  >  0, 
preprocess  S  so  that  for  any  query  triangle,  we  can  quickly  compute  the  number  of 
points  contained  in  that  triangle.  We  give  an  algorithm  with  0(--7i^log^''^n)  query 
time,  using  0{m)  space.  The  preprocessing  time  is  bounded  by  0(n>/m  log'^"'"^/^  n). 
Similar  bounds  have  been  obtained  independently  by  Chazelle  [Che]. 

(x)  Overlapping  planar  maps:  Given  two  planar  maps  P,  Q,  and  a  bivariate  function 
Fp{x,y),  Fg(x,y)  associated  with  each  of  them,  such  that  over  each  face  of  P  the 
function  Fp  has  some  simple  structure  (e.g.  it  is  constant,  linear,  or  convex  over  each 
face),  and  similarly  for  Q,  determine  a  point  that  minimizes  Fp{x,y)  -  FQ{x,y). 
We  show  that  if  the  maps  satisfy  certain  conditions,  then  an  optimal  point  can  be 
computed  in  0{n*f^  log^'^+^^Z^  n)  time,  where  n  is  the  total  complexity  of  the  two 
maps.  The  details  of  this  application  can  be  found  in  [Age]. 
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Our  algorithm  for  partitioning  the  plane  works  in  two  phases.  The  first  phase  partitions 
the  plane  into  0(r^  log""  r)  triangles  so  that  no  triangle  intersects  more  than  O(^)  hnes. 
The  second  phase  reduces  the  number  of  triangles  to  0{r'^)  till  maintaining  the  property 
that  no  triangle  meets  more  than  O(^)  lines,  using  a  similar  technique  to  that  of  Matousek 
[Maa].  We  beheve  that  the  second  phase  is  redundant,  i.e.  that  the  first  phase  or  some 
appropriate  variant  of  it  produces  only  O(r^)  triangles,  but  we  have  not  been  able  to  prove 
it  so  fax.  Even  if  this  is  not  the  case,  in  most  of  the  applications  we  can  stop  after  the 
first  phase  and  solve  the  subproblems  directly  within  each  of  the  resulting  0{r-^  log"^  r) 
triangles,  without  increasing  the  asymptotic  complexity  of  the  algorithm. 

The  paper  is  organized  as  follows.  In  Section  2  we  describe  the  geometric  concepts  that 
we  use  in  this  paper.  Section  3  gives  an  algorithm  for  a  subproblem  which  is  interesting  in 
its  own  right,  namely  that  of  computing  the  k"*  leftmost  intersection  point  induced  by  a 
set  of  n  lines  inside  a  given  convex  quadrilateral.  Section  4  is  the  heart  of  our  algorithm; 
it  presents  an  algorithm  to  partition  a  given  convex  quadrilateral  Q  intersecting  m  lines 
and  containing  K  of  their  intersection  points  into  0{j-  +  -t^)  convex  quadrilaterals  (for  an 
arbitrary  parameter  C)  so  that  both  the  nimiber  of  Unes  crossing  any  subquadrilateral  Q' 
and  the  number  of  their  intersections  inside  Q'  are  small  (the  precise  conditions  are  given 
in  Theorem  4.15).  In  Section  5  we  describe  the  first  phase  of  the  partitioning  algorithm, 
which  basically  consists  of  applying  the  algorithm  of  Section  4  recursively,  and  Section  6 
describes  the  second  phase  of  the  algorithm.  Section  7  shows  how  to  modify  the  analysis 
to  handle  degenerate  cases  as  well.  In  Section  8  we  conclude  with  some  final  remarks  and 
open  problems. 

2      Geometric  preliminaries 

This  section  defines  the  geometric  concepts,  and  formalizes  the  notation  that  we  will  be 
using  in  this  paper.  Let  H  =  {^1,^2,  •  •  •  ,^n}  denote  a  set  of  n  lines  in  the  plane.  These 
lines  induce  a  planax  map  called  the  arrangement  A{H)  of  H,  whose  vertices  axe  the 
intersection  points  of  lines  in  Ti,  edges  axe  maximal  connected  portions  of  Unes  in  Ti  not 
containing  a  vertex,  suid  faces  axe  maximal  connected  portions  of  the  plane  not  meeting 
any  edge  or  vertex  of  ACH).  See  [Ed]  for  more  details.  For  simplicity  we  assume  that  the 
Unes  in  7{  axe  in  general  position,  that  is,  no  three  Unes  axe  concurrent,  and  that  no  line 
in  7i  is  vertical.  For  any  point  p  €  R^,  we  define  its  level  to  be  the  number  of  Unes  in  H 
lying  above  it  (not  counting  the  Unes  passing  through  p).  For  any  0  <  k  <  n,  the  k-level 
of  A{7i)  is  the  set  of  edges  of  ACH)  whose  level  is  k  (see  figure  1).  A  A;-level  of  ACH)  is 
a  monotone  polygonal  chain  with  two  unbounded  rays.  More  details  on  Jfc-levels  can  be 
found  in  [Ed],  [EW],  [Wea]. 

Lemma  2.1  Letp\  andp2  be  two  points  in  the  plane  whose  levels  are  ki  and  ^2  respectively. 
Then  the  line  segment  pipi  intersects  at  least  \ki  —  k2\  lines  ofH. 
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plified  2-level 

^1-approx  2-level 
?^2-level 


Figure  1:  2-level,  1-approximate  2-level,  2-simplified  2-level 


Proof:  Obvious  from  the  definition  of  fc-levels. 


a 


Lemma  2.2  ([Maa])   Let  U\M2i  ■■■  Ms  denote  ■pair-wise  disjoint  seta  of  levels  of  A{C) 


n 


with  \Ui\  >  u!  for  all  i  <  s,  then  there  exists  a  level  U,  €  Ui  such  that  J2i=i  \U,\  <  —■ 


We  call  an  x-monotone  polygonal  path  11  (not  necessarily  formed  by  the  edges  of  ACH)) 
an  €- approximate  k-level,  for  e  <  A;,  if  it  lies  in  the  strip  bounded  between  the  k  —  e  and 
k  +  e  levels  of  AiH)  (see  figure  1).  A  set  of  e-approximate  2ez-levels,  for  i  <  [^J ,  is  called 
an  e-approximate  leveling  of  ACH).  Let  po,p\,  ...  ,Pm  be  the  vertices  of  a  A;-level  of  ACH). 
For  any  6  <  m,  we  define  the  8-simplified  k-level  to  be  the  polygonal  path  formed  by 
connecting  po  to  p^,  ps  to  p^s,  ■  ■ .,  P[^j5  to  pm,  and  concatenated  with  the  left  and  right 
rays  of  the  fc-level  incident  to  p\  and  pm  respectively.  Edelsbrunner  and  Welzl  [EW]  proved 
that 

Lemma  2.3  ([EW])   For  a  given  set  of  n  lines  H.,  and  0  <  fc  <  n,  a  8-simplified  k-level 
is  a  \-^- approximate  k-level  of  A^H). 

U 

Let  Q  be  a  convex  quadrilateral  with  vertices  91,921 93?  94  ordered  in  counter-clockwise 
direction,  and  let  dQ  denote  the  boundary  of  Q.  Let  £  =  {^1,^2,  •  •  •  ,^m}  denote  a  set  of 
lines  passing  through  Q,  and  let  K  be  the  number  of  intersection  points  of  C  contained 
inside  Q.  The  intersection  points  of  dQ  and  C  are  called  the  endpoints  of  C.  For  simplicity 
let  us  assume  that  every  fine  of  £  intersects  dQ  at  two  points,  i.e.  does  not  touch  Q  at 
one  of  its  vertices  (or  at  an  edge). 

The  set  £  of  lines  can  be  partitioned  into  three  different  classes  according  to  the  location 
of  their  endpoints  (see  figure  2). 
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Figure  2:  Q  and  C:  solid  lines  denote  red  lines,  dashed  lines  denote  green  lines  and  dotted  lines 
denote  blue  lines. 

(i)  "Red"  lines  Cr-  a  line  i  belongs  to  Cr  if  one  of  its  endpoints  lies  on  q^  or  q^ 
(referred  to  as  the  left  endpoint)  and  the  other  one  lies  on  gTgl  or  ^3^  (referred  to 
as  the  right  endpoint),  i.e.  £  crosses  the  diagonal  5T53;  let  |£r|  =  mr- 

(ii)  "Blue"  lines  £(,:  a  Une  i  belongs  to  Cb  if  one  of  its  endpoints  lies  on  q^  (referred  to 
as  the  left  endpoint)  ajid  the  other  one  lies  on  92^  (referred  to  as  the  right  endpoint); 
let  \Cb\  =  TTib- 

(iii)  "Green"  lines  Cgi  a  line  £  belongs  to  Cg  if  one  of  its  endpoints  lies  on  qiq^  (referred  to 
as  the  left  endpoint)  and  the  other  one  lies  on  gj^  (referred  to  as  the  right  endpoint); 
let  \Cg\  =  Trig. 

These  three  sets  axe  pairwise  disjoint  and  Cr^CgUCb  =  C.  Let  K^y  denote  the  number 
of  intersection  points  between  Cx  and  Cy  lying  inside  Q,  where  x,y  €  {r,b,g}.  It  is  easily 
seen  that  Cb  and  Cg  do  not  intersect  inside  Q,  that  is  Kbg  =  0.  Without  loss  of  generality 
we  can  assume  that  rrir  >  y  because  otherwise  we  can  cyclically  renumber  the  vertices  of 
Q  so  that  the  above  inequality  is  satisfied.  For  the  saJke  of  exposition  we  assume  that  51  is 
the  upper  left  comer  of  Q. 

Let  A  =  {oi,  ...  ,am,}  (resp.  B  —  {61,  ...  ,6mr})  be  the  left  (resp.  right)  endpoints 
of  Ct  appearing  in  counter-clockwise  (resp.  clockwise)  direction  around  dQ.  For  each  line 
£i  G  Cr  let  a^(j)  (resp.  &<T(t))  denote  its  left  (resp.  right)  endpoint. 

3      Selecting  the  «:*^  leftmost  intersection  point 

Let  pi ,  ...  ,pk  denote  the  intersection  points  of  the  Hnes  in  C  lying  inside  Q  and  ordered 
in  increasing  i-coordinate.  In  this  section  we  consider  the  problem  of  calculating  the  k"* 
leftmost  point  p^  of  these  intersections,  for  any  given  k  <  K.  The  algorithm  extends  the 
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recent  algorithm  of  Cole  et  al.  [CSSS]  for  calculating  the  A:"*  leftmost  intersection  point  of 
n  lines  in  the  entire  plane.  One  of  the  key  steps  of  the  algorithm  is  to  count  the  number  of 
intersection  points  of  C  lying  inside  a  given  convex  quadrilateral.  We  first  give  an  algorithm 
that  counts  the  number  of  intersections  inside  a  convex  quadrilateral  in  time  0(m  log  m), 
where  m  is  the  number  of  lines  passing  through  Q. 

3.1      Counting  intersection  points  inside  a  convex  quadrilateral 

We  partition  the  set  of  Hnes  £  into  £,,  £&  and  Cg  as  defined  earlier,  and  count  each  of 
Kxy,  for  x,y  G  {r,b,g},  separately.  First,  we  show  how  to  count  the  number  of  red-red 
intersection  points  lying  inside  Q,  i.e.  Krr-  It  is  obvious  that  two  red  lines  i,  and  £j  in  C 
intersect  if  and  only  if  7r(z)  <  7r(;)  and  a{i)  >  (x{j),  or  vice  versa.  We  thus  obtain  Krr 
by  counting  the  number  of  inversions  in  u  with  respect  to  tt,  which  is  easily  done  in  time 
0{mr  logmr)  by  sorting  the  endpoints  of  Cr  adong  the  boundary  of  Q.  Similarly,  we  can 
compute  Kbb  (resp.  Kgg)  in  time  0(m6  log  mfc)  (resp.  0{mg\ogmg)). 


Figure  3:  Green  line  £  and  its  endpoints  61,62;  ^1,  £4  do  not  intersect  i 


As  for  determining  Krg,  let  ei  (resp.  62)  be  the  left  (resp.  right)  endpoint  of  a  green 
line  £.  The  Une  £  intersects  a  red  line  £j  if  the  right  endpoint  of  ij  lies  on  either  e^  or  q^ 
(see  figure  3).  Therefore  using  the  same  sorted  Ust  B,  we  can  compute  the  number  of  red 
lines  that  £  intersects  by  locating  Ci  and  62  in  B,  which  impUes  that  Krg  can  be  computed 
in  time  0{mglogmr).  Similarly,  Krb  can  be  computed  in  time  0{mblogmr).  Thus,  we 
have  shown 

Lemma  3.1   For  a  given  convex  quadrilateral  Q  and  a  set  C  of  m  lines  passing  through 
it,  we  can  count  the  number  of  intersections  of  C  lying  inside  Q  in  time  O(mlogm). 
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Remark  3.2:  The  above  algorithm  can  be  extended  to  count  the  number  of  intersections 
lying  inside  any  Jt-gon  B  in  0(km\ogTn)  time,  by  triangulating  B,  and  by  counting  the 
number  of  intersections  in  each  triangle  separately. 

Next  we  give  an  algorithm  for  finding  the  «"*  leftmost  intersection  point  inside  Q  using 
the  counting  procedure  just  described. 

3.2      Computing  the  k*^  leftmost  intersection 

Cole  et  al.  [CSSS]  have  given  an  0(m  log m)  algorithm  that  returns  the  /c*'*  leftmost  inter- 
section point  in  a  set  £  of  m  lines  in  the  plane.  Their  algorithm  cannot  be  applied  directly 
to  our  case  because  it  counts  all  intersection  points  of  C  lying  to  the  left  of  the  k^^  one, 
while  we  count  only  those  intersection  points  that  lie  inside  Q.  However,  we  will  show  that 
the  algorithm  of  [CSSS]  can  be  eaisily  adapted  to  our  case.  We  assume  that  the  reader  is 
familiau:  with  this  algorithm,  so  we  describe  it  here  only  briefly,  and  show  how  to  modify 
it  for  our  case;  for  more  details,  see  [CSSS]. 

Let  7r(a)  denote  the  permutation  of  hnes  in  £  sorted  in  decreasing  order  of  their  inter- 
cepts with  the  Hne  x  =  a,  assuming  that  no  two  lines  in  £  intersect  at  i  =  a.  The  rank 
of  a,  denoted  by  ?/'(a),  is  defined  to  be  the  number  of  mtersection  points  of  £  lying  inside 
Q  and  to  the  left  of  i  =  a.  The  k"*  leftmost  intersection  point  p^  is  such  that  i/'(a*)  =  k, 
where  a'  =  x{p^)  +  e,  for  some  sufficiently  small  e  >  0  (in  this  description  we  assume  that 
no  two  points  p,,  pj  have  the  same  z-coordinate). 

Cole  et  aJ.'s  algorithm  [CSSS]  is  based  on  Megiddo's  technique  [Me],  further  improved 
by  Cole  [Co],  for  implicit  parametric  searching.  Using  a  sorting  network,  e.g.  that  of  Ajtai 
et  al.  [AKS],  7r(a*)  can  be  obtained  by  performing  0(m log m)  comparisons,  performed  in 
O(log  m)  batches  of  0{m)  independent  comparisons  each,  where  each  comparison  answers 
a  question  of  the  form:  "Does  ii  he  above  £j  at  x  =  a*?".  Although  we  do  not  know  a",  we 
can  answer  a  question  of  this  form  by  computing  the  rank  of  the  abscissa  x{qij),  where  g,j 
is  the  intersection  point  of  £i  and  £j.  If  %l;(x{qij))  <  k,  then  g,j  Hes  to  the  left  of  x  =  a*,  so 
the  order  of  £{  and  ij  there  coincides  with  their  order  at  x  =  — oo;  otherwise  the  order  is 
reversed.  In  such  a  sorting  network  there  are  O(logm)  levels  and  at  each  level  s  =  0{m) 
questions  axe  asked;  let  Zi,Z2,  ...  ,2,  denote  the  x-coordinates  of  the  intersection  points 
corresponding  to  the  questions  to  be  answered  at  a  level  j.  If  we  resolve  the  question  for 
the  median  value  Zmed  oi  zi,  ...  ,2,,  then  we  can  resolve  the  question  for  at  least  half  of 
the  above  points,  e.g.  if  z^ed  <  a',  then  Zi  <  a'  for  aW  Zi  <  z^ed-  For  the  remaining  half 
of  the  points,  the  question  can  be  resolved  by  applying  the  same  method  recursively,  thus 
allowing  us  to  resolve  all  m  questions  by  actually  computing  the  rank  of  only  (9(log  m) 
points.  Once  all  comparisons  are  resolved  at  level  j,  we  can  continue  to  level  j  +  1  oi  the 
network.  At  the  end,  all  comparisons  have  been  resolved  and  we  have  obtained  the  sorted 
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order  of  the  lines  in  £  along  x  =  a',  still  without  knowing  the  exact  value  of  a'.  However, 
this  order  is  easily  seen  to  determine  uniquely  how  many  intersection  points  of  C  lie  to 
the  left  of  I  =  a',  thus  the  set  of  these  intersections,  as  well  as  the  subset  of  these  points 
lying  inside  Q,  are  fully  determined  at  the  end  of  the  algorithm.  By  keeping  track  of  all 
inequalities  of  the  form  z^^d  <  a'  or  z^ed  >  a'  we  finally  obtain  an  interval  a  <  a'  <  /3, 
and  the  logic  of  the  procedure  is  easily  seen  to  imply  that  a  is  the  x-coordinate  of  p«.  The 
point  Pk  itself  is  also  easy  to  obtain  by  recording,  for  each  ^med,  the  two  lines  that  induce 
it. 

To  resolve  a  comparison  at  x  =  ^med,  we  proceed  as  follows.  Let  Q~  denote  the  portion 
of  Q  lying  to  the  left  of  the  vertical  hne  x  =  Zmed-  by  Lemma  3.1  and  the  subsequent 
Remaxk,  we  can  count  the  number  of  intersection  points  of  C  lying  inside  Q~  (which  is 
at  most  a  pentagon)  in  0(m  log  tti)  time.  Thus  the  rank  of  ^med  can  be  computed  in 
0(m log  m)  time,  so  the  overall  running  time  of  the  algorithm  is  0(mlog~'m).  By  using 
Cole's  improvement  [Co],  the  time  can  be  reduced  to  O(mlog^m).  Cole  et  al.  [CSSS] 
observed  that  it  is  not  necessaxy  to  compute  the  exact  rank  of  a  point  because  we  are  only 
interested  in  knowing  whether  Zj^td  lies  to  the  left  or  to  the  right  of  p«.  They  proved  that  it 
suffices  to  compute  an  approximate  rank  of  a  point,  which  can  be  done  in  0{m)  amortized 
time,  yielding  an  0(m  log  m)  algorithm.  Using  the  same  idea  we  can  also  compute  our 
version  of  approximate  rank  in  amortized  time  0{m),  and  therefore  we  obtain 

Theorem  3.3  Given  a  convex  quadrilateral  Q,  a  set  of  m  lines  C  and  an  integer  k,  we  can 
determine,  in  time  0{mlogm),  a  vertical  line,  so  that  Q  contains  exactly  k  intersection 
points  to  the  left  of  that  line. 

D 

Remark  3.4: 

(i)  The  remark  following  Lemma  3.1  impUes  that  we  can  extend  the  above  algorithm  to 
any  k-gon  with  0(A;m  log  m)  rvmning  time. 

(ii)  The  above  algorithm  uses  the  sorting  network  of  Ajtai  et  al.  [AKS],  which  involves  a 
very  large  constant,  therefore  making  our  algorithm  impractical.  A  possible  solution 
to  circumvent  this  problem  is  to  use  the  much  simpler  Batcher's  sorting  network 
[Ba];  this  network  has  depth  O(log^  m)  and  therefore  the  overall  running  time  of  the 
resulting  algorithm  becomes  O(mlog^m).  However,  to  obtain  a  better  asymptotic 
complexity,  we  will  be  using  the  former  algorithm. 

4     Partitioning  a  convex  quadrilateral 

This  section  describes  an  algorithm  to  partition  the  convex  quadrilateral  Q  into  convex 
quadrilateral  subcells  so  that  on  the  average  only  few  lines  pass  through  each  cell,  and  the 
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maximum  number  of  lines  passing  through  each  cell  is  also  low.  To  be  more  precise,  for 
any  given  positive  integer  (  <  m,  we  want  to  partition  Q  into  M  =  0{^  +  ■^)  cells  so  that 
the  following  conditions  are  satisfied  (where  m,  is  the  number  of  lines  passing  through  the 
z"*  cell,  and  Ki  is  the  number  of  intersections  within  the  i"*  cell): 

(i)      ^m,  =  O  [m  +  —  J  ,      maxm,  =  O(V^) 

'm 
(ii)     ^i^.  =  iiT,  max/sT,  <  C^K. 


4.1      A  special  case 

Before  solving  this  general  problem  we  consider  a  special  case  in  which  Q  contains  only 
"red"  lines  and  we  do  not  care  how  many  intersection  points  lie  in  a  cell.  Moreover, 
although  most  cells  in  the  construction  to  follow  will  be  quadrilaterals,  some  could  have 
up  to  six  edges. 


Figure  4:   Partitioning  Q  into  Qi,  ...  ,Qt,  with  m  =  6  and  C  =  2.  The  dashed  edges  are  pseudo 
edges. 

First,  compute  the  left  and  right  endpoints  of  all  lines  in  £  and  sort  them  to  obtain 
the  Usts  A  and  B.  Next,  for  all  1  <  i  <  <  =  \^]  connect  a,-^  to  bi(.  The  segments  a^^  are 
called  pseudo  edges  (see  figure  4).  The  pseudo  edges  partition  Q  into  cells  Qi,  Q2,  ...  ,  Q«. 
The  following  sequence  of  lemmas  show  that  this  partitioning  has  the  desired  properties. 

Lemma  4.1  A  line  ii  6  C  intersects  at  least  \n(i)  —  a{i)\  lines  of  C  inside  Q. 
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Proof:  Without  loss  of  generality  assume  that  7r(t)  >  a{i).  Since  only  a(i)  -  1  lines  have 
their  right  endpoints  preceding  the  right  endpoint  of  £,,  at  least  7r(z)  -  a{i)  of  the  hnes, 
whose  left  endpoints  precede  that  of  £,,  have  their  right  endpoints  after  b„(i)  in  B.  But 
every  such  hne  crosses  the  line  i,  inside  Q,  showing  that  at  least  \n{i)  -  a{i)\  Unes  of  £ 
intersect  ii  inside  Q. 

a 

Let  6i  denote  the  number  of  pseudo  edges  intersected  by  the  Hne  £,. 

"*  2K 

Lemma  4.2  ^  ^i  <  — : — 1-  "i- 

t=i  ^ 

Proof:  If  a  hne  ii  intersects  6,  pseudo  edges,  then  it  is  easily  seen  that 

|7r(0  -  cr{i)\   >  (S,  -  1)C. 
Summing  over  all  lines,  we  obtain 

m  m 

i=l  i=l 

Let  Vi  denote  the  nxunber  of  lines  of  £  that  £,  intersects  inside  Q.  We  have  YlT=i  ^>  = 
2i\,  because  there  are  only  K  intersection  points  inside  Q  and  each  intersection  is  counted 
twice.  By  Lemma  4.1,  |7r(i)  —  cr{i)\  <  Ui,  therefore 

J2{6,-1)C<2K     or     5^6,  <  —  +  m. 

i=l  i=l  ^ 

D 

Let  m,  denote  the  number  of  lines  passing  through  the  cell  Q,.  Since  the  line  £.j  lies  in 
6j  +  1  cells,  we  have 

t  m  /TV-  \ 

Corollary  4.3  J2^i  =  J2i^J  +  1)<2—  +  m  □ 

t=i  j=i  V  C        / 

Lemma  4.4  Let  fii  denote  the  number  of  lines  passing  through  the  cell  Q,  but  not  having 
any  of  their  endpoints  on  dQi.   Then  fii  <  2VK,  for  all  i  <t. 

Proof:  Let  £,  denote  the  set  of  Unes  passing  through  the  cell  Q,  but  not  having  any 
endpoint  in  Qi.  Partition  the  set  £,  into  disjoint  subsets  C^y,  for  x  =  I,  . . .  ,t  —  i,  y  = 
1,  ...,?  —  1,  such  that  a  line  ij  belongs  to  Cxy  if  one  of  its  endpoint  lies  in  Qi+x  and  the 
other  in  Qi-y.  Let  fi^y  =  \£xy\-  Note  that 

V^,  €  £xv,  K(j)  -  a{j)\   >   C(x  +  y  -  1). 
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Figure  5:  ij  intersects  at  least  6C  lines  of  C. 


By  Lemma  4.1  a  line  ij   G  Cry  intersects  at  least  (x  +  y  —  1)C  lines  (see  figure  5). 
Summing  over  all  i  and  y  we  get 


2^  ^  EE/^-v(^  +  y-i)-c 

2K 


X      y 


or     —     >     X^l^a:-/i^y  +  ^^y /x^y  -  ^^^xy 


a;      y 


Let  C^  C  £j  denote  the  set  of  all  lines  having  one  endpoint  in  Qi+x-  Then 

X      y  X  y  X 

Similarly,  if  Cy  C  £j  denotes  the  set  of  hnes  having  one  endpoint  in  Q,_y,  then 

H^y /^xy  =  Y>y\^y\- 

X      y  y 

Moreover  J2x  Yly  y-xy  =  A*t-  Therefore  (4.1)  can  be  written  as 


E^-i^l  +  Lyl-^yl-M.  < 


2K 


X  y  S> 

Since  $2j.  |£iI  =  Mt  3^<i  for  all  X,  |£!^|  <  2C,  the  term  Hx  ^l^il  minimizes  when 


141  = 


2C 

/i,-2C 
I    0 


L2C, 


for  1  <  X  < 

forx=^- 
otherwise 


2C 
+  1 


(4.1) 


(4.2) 


(4.3) 
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Hence, 


LffJ 


E^-l^xl  >   E^-2C+  U-2C 


r=l 


tJ-x 


2C 


L2C 


>  HC 

-      2 


2CJ 


+  1 


.2C 


Similarly, 


Q|J,,),(„_,C[|J)([|J..) 
-0(''-4|J) 


2C 


2CJ 


y 
By  substituting  (4.4)  and  (4.5)  in  (4.3),  we  obtain 


Let 


2C 


=  a  and  ^  =  a  +  i/,  for  some  0  <  i/  <  1,  then  (4.6)  becomes 

IK 

^      >     2(a  +  1)  (2C(a  +  j/)  -  (a)  -  2C(a  +  v) 

K 


or     —      >      (a  +  l)(a  +  2i/)  -  (a  +  i/) 


C^ 


=     a^  +  2ai/  +  V 

=     (a  +  i/)^  +  1/(1  -  v) 

>     (a  +  i/)^  (because  0  <  i/  <  1) 

ill 
4^2 


Therefore,      /x^  <  2\/^. 


[4.4) 


(4.5) 


(4.6) 


Concerning  running  time,  notice  that,  the  sets  A  and  B,  and  therefore  Qi,  ...  ,  Q(  can 
be  computed  in  time  0{m  log  m).  Moreover,  for  each  Une  i.j  we  can  easily  compute  the  cells 
it  crosses —  they  axe  simply  all  the  cells  lying  between  (and  including)  the  cell  containing 
a^(j)  and  the  cell  containing  b„(j) .  Thus  the  total  time  spent  in  computing  the  lines  passing 
through  each  cell  is  bounded  by  0{m log m  +  K/C,).  Now  it  follows  from  Theorem  3.3  that 
number  Ki  of  intersections  within  Q,  can  be  computed  in  time  0(mi  log  m,).  Therefore 
the  total  time  spent  in  computing  intersections  is  at  most 

OiY^rriilogmA    =  of  (^m.)logmj    =   Onm  +  — )logmj. 
Hence,  we  can  conclude  that 
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Lemma  4.5  The  above  algorithm  partitions  Q  into  t  =  \y]  subcells  Qi,  ■  ■ .  ,Qt  in  time 
0(m  log  m)  such  that  (i)  Yll^i  m,  <2[m-\-  j^,  and  (ii)  max,  m,  <  2\/F  +  2C.  Moreover, 
■we  can  compute  Ci  and  Ki  (as  defined  above)  for  all  Q,,  in  time  0{{m  +  — )logm). 

D 

Remark  4.6:  Note  that  in  the  special  case  at  hand  we  have  ignored  the  issue  of  making 
the  numbers  K,  small.  Also,  the  resulting  partitioning  can  have  one  cell  with  6  edges 
or  two  cells  with  5  edges  each.  These  issues  will  be  addressed  in  the  following  general 
algorithm. 

4.2      General  Algorithm 

Next  we  give  an  algorithm  for  the  general  case  using  the  above  procedure  as  a  subroutine. 
Our  algorithm  consists  of  five  steps. 


Figtire  6:  Red,  green  and  blue  pseudo  edges 

The  first  step  applies  the  preceding  algorithm  to  the  collection  of  red  Unes,  to  partition 
the  quadrilateral  Q  into  a  collection  of  cells  Q""  =  {Q\,  ■  ■ .  ,Ql},  for  t  =  [^] ,  so  that  every 
cell  contains  the  left  (resp.  right)  endpoints  of  at  most  C  red  hnes.  Let  R  =  {ri,  ...  ,  rt-i} 
denote  the  set  of  pseudo  edges  that  bound  the  quadrilaterals  of  Q*";  they  will  be  referred 
to  as  "red"  pseudo  edges. 

The  second  step  appUes  the  preceding  algorithm  to  the  collection  of  green  lines  (where 
this  time  dQ  is  divided  into  a  left  and  a  right  portions  at  54  amd  ^2),  to  partition  Q  into 
a  collection  of  cells  Q^  =  {Q{,  . . .  ,Ql},  for  u  =  f^] ,  so  that  every  cell  contains  the  left 
(resp.  right)  endpoints  of  at  most  C  green  lines.  Let  G  =  {firi,  ...  ,gu-i}  denote  the  set 
of  pseudo  edges  bounding  the  quadrilaterals  of  Q^;  these  pseudo  edges  will  be  referred  to 
as  "green"  pseudo  edges.  Since  the  endpoints  of  all  green  Unes  lie  on  q^  and  5354,  all 
green  pseudo  edges  extend  from  qiq4  to  9394.  Assume  that  the  quadrilaterals  Qf  are  sorted 
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in  their  order  along  9^4,  and  that  Qf  is  bounded  by  the  two  pseudo  edges  y,_i,  g,,  for 
1  <  i  <  u. 

In  the  third  step  we  apply  the  same  algorithm  to  the  collection  of  blue  lines  (again, 
partitioning  dQ  at  92  and  54),  to  partition  Q  into  a  collection  of  cells  Q*  =  {Q\,  ...  ,  Q*}, 
for  V  =  [^^1,  so  that  every  cell  contains  the  left  (resp.  right)  endpoints  of  at  most  (  blue 
lines.  Let  B  =  {bi,  . . .  ,6v-i}  denote  the  set  of  pseudo  edges  bounding  the  quadrilaterals 
of  Q*;  they  are  referred  to  as  "blue"  pseudo  edges.  This  time  all  pseudo  edges  extend  from 
q^  to  q^.  Assume  that  the  quadrilatercds  Q-  are  sorted  in  their  order  along  qiq2,  and 
that  Q-  is  bounded  by  the  two  pseudo  edges  6i_i,  6i,  for  1  <  i  <  u.  See  figure  6  for  an 
illustration  of  the  resulting  pseudo  edges. 


Figure  7:  Internal  and  boundary  cells:  shaded  regions  denote  internal  cells 


Let  S  =  R\J  B  U  G,  and  S  the  set  of  endpoints  of  segments  in  S.  Let  F  denote  the  set 
of  bounded  faces  in  the  arrangement  of  dQ  U  S.  If  a  face  ^q  £  F  touches  the  boundary 
of  Q,  then  it  is  called  a  boundary  cell,  otherwise  it  is  called  an  internal  cell  (see  figure  7). 
The  boundary  of  an  internal  cell  is  formed  by  four  pseudo  edges,  two  of  which  axe  aJways 
red.  On  the  basis  of  the  color  of  their  edges  we  partition  the  cells  into  three  categories: 

(i)  red-green  cells:  cells  formed  by  the  intersection  of  Qf  for  1  <  i  <  u  and  Q^  for 
!<;<*•  A  red-green  cell  formed  by  Qf  PI  Q^  is  denoted  by  :Fij.  If  Tij  is  an  internal 
cell,  then  its  edges  are  portions  of  gi_i,  g,  G  G  and  r_,_i,  rj  6  R  (and  i  <  u,  1  <  j  <  t). 

(ii)  red-blue  cells:  cells  formed  by  the  intersection  of  Q-  for  1  <  i  <  u  and  Q^  for 
1  <  j  <  <•  A  red-blue  cell  formed  by  Q-  fl  Q^  is  denoted  by  Tij.  If  Tij  is  an  internal 
cell,  then  its  edges  are  portions  of  6j_i,  bi  ^  B  and  r_,_i,  Vj  £  R  (and  i  <  v,l  <  j  <  t). 
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(iii)  red- blue- green  cells:  cells  formed  by  the  intersection  of  Q{,  Qj  and  Q-  for  1  <  j  <  t. 
We  denote  such  a  cell  by  !Fu-  If  ^u  is  an  internal  cell,  then  its  edges  are  portions  of 
bi,  gi,  r,_i  and  r,  (and  1  <  i  <  t). 

The  problem  with  the  boundary  cells  produced  by  the  arrajigement  of  dQ  U  <5  is  that 
some  of  them  may  have  more  than  four  edges.  However,  since  every  cell  has  at  most  four 
pseudo  edges,  it  cannot  have  more  than  eight  edges  in  total.  At  the  fourth  step  of  our 
algorithm,  we  partition  all  boundary  cells  with  more  than  four  edges  into  two  or  three 
convex  quadrilaterals  and  triangles  by  appropriate  diagonals,  which  we  also  call  pseudo 
edges  (see  figure  8). 


^^^^ 

^ 

y^           ^ 

■v 

-»-^^l 

■  V    ■' 

n- 

-1  .  ■     .■•' 

■<7fl 

•i 
> 

/ 

b: 

\ 

■N, 

-s/ 

/r. 

Figure  8:  Cells  with  seven  or  eight  edges:  djished  lines  are  the  pseudo  edges  added  to  partition 
the  cell  into  convex  quadrilaterals  and  triangles 

A  cell  with  five  or  six  edges  is  partitioned  arbitrarily  into  two  convex  quadrilaterals, 
but  a  cell  with  seven  or  eight  edges  is  partitioned  in  such  a  way  so  that  one  of  the  new 
pseudo  edges  intersects  only  red  and  green  Unes,  and  the  other  intersects  only  red  and  blue 
hnes.  It  is  easily  seen  that  such  a  partition  is  always  possible.  For  example,  an  8-cell  lying 
between  the  red  pseudo  edges  r,_i,  tj  is  partitioned  by  connecting  the  left  endpoint  of  r,_i 
to  the  left  endpoint  of  r,,  and  the  right  endpoint  of  r,_i  to  the  right  endpoint  of  r,  (see 
figure  8).  7-cells  axe  partitioned  in  a  similar  manner. 

Note  that  there  can  be  at  most  one  cell  having  7  or  8  edges.  The  existence  of  an  8-cell 
implies  that  no  pair  of  pseudo  edges  intersect,  and  the  existence  of  a  7-cell  allows  only 
very  limited  pattern  of  pseudo-edge  intersections.  Moreover  the  existence  of  a  7  or  8-cell 
impUes  that  the  first  step  creates  at  least  one  red  pseudo  edge. 

Next,  for  each  of  the  cells  produced  so  far,  we  compute  the  lines  passing  through  it, 
and  count  the  number  of  intersection  points  lying  in  it.  If  i^  >  6C^  and  a  cell  Ta  has  Ka 
>  C,\/K  intersection  points,  then  we  partition  it  into  |"^^]  cells,  each  containing  at  most 

C,\/K  intersections,  using  the  algorithm  of  Section  3.  With  some  caxe  we  can  obtain  such 
a  partition  with  at  most  one  5-cell,  which  we  then  split  by  a  diagoned  into  two  subcells. 
This  yields  a  final  collection  of  quadrilaterals  [Qi,  Q2,  ...  ,  Qm},  which  is  the  output  of 
the  algorithm. 
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4.3      Analysis  of  the  algorithm 

The  algorithm  just  described  for  partitioning  a  convex  quadrilateral  is  fairly  simple.  Its 
analysis,  however,  is  not.  For  i  =  1,  ...  ,  M,  let  m,  denote  the  number  of  Hnes  passing 
through  Qi,  and  let  K,  denote  the  number  of  intersections  contained  in  Q,.  For  simplicity 
we  assume  that  the  lines  are  in  general  position,  that  is,  no  three  hnes  axe  concurrent,  and 
no  line  is  vertical.  We  begin  by  bounding  the  total  number  of  cells  created  by  the  algorithm 
and  the  number  of  lines  meeting  each  cell.  As  we  shall  see  later,  the  constants  appearing 
in  the  bounds  for  max,  m,-  and  ^,  m,  control  the  exponent  of  the  logarithmic  factor  in  the 
bound  for  the  time  complexity  of  the  overall  paxtitioning  algorithm,  therefore  we  try  to 
obtain  as  small  constants  as  possible.  If  we  do  not  worry  about  constants,  the  proofs  can 
be  simplified  a  lot.  To  bound  the  total  nim[iber  of  cells  produced  by  our  algorithm,  we 
first  estimate  the  number  of  cells  formed  by  overlapping  the  original  red,  green  and  blue 
pseudo  edges. 


Lemnia  4.7  Let  N  be  the  number  of  cells  produced  by  overlapping  the  red,  green  and  blue 
pseudo  edges,  then 

N  <   ^m  +  ^  -  2.  (4.7) 


^2  K       ^ 

<   ^m  +  -  -  2. 


Proof:  We  bound  the  number  of  boundary  and  internal  cells  separately.  The  number  of 
boundary  cells  is  obviously  equal  to  the  number  of  endpoints  of  pseudo  edges,  namely  2|<5| 
=  2{u  +  V  +  t  -  3). 
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Figure  9:  Charging  of  a  red-green  intersection  point 


Next,  consider  the  number  of  internal  cells.  To  bound  the  number  of  red-green  and 
red- blue- green  cells,  we  use  the  following  charging  scheme.  Let  w  be  the  intersection  point 
of  a  red  hne  t  and  a  green  Hne  i' .  Suppose  that  the  right  endpoint  of  i  hes  on  ^T^J  in  a 
cell  QJ,  and  that  the  left  endpoint  of  £'  (on  gT^)  hes  in  Q,\.  If  i  =  1  and  1  <  j  <  «,  then 
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we  charge  w  to  the  red-blue-green  cell  J^ij  (e.g.  W4  in  figure  9)  ,  and  if  1  <  i  <  u  and 
I  <  j  <  t,  then  we  charge  w  to  the  red-green  cell  ^,j  (e.g.  W2  in  figure  9).  If  Qf  n  Q^  is 
a  boundary  cell,  then  we  do  not  charge  w  to  any  cell  (e.g.  W3  in  figure  9).  Similarly,  if 
the  right  endpoint  of  £  Hes  on  gigj,  then  we  charge  w  to  the  cell  J^,j  (e.g.  lOj  in  figure  9), 
where  Qf  contains  the  right  endpoint  of  £'  (on  53^4)  and  Q^  contains,  as  before,  the  right 
endpoint  of  £.  Again,  no  charge  is  made  if  Qf  fl  Q^  is  a  boundary  cell. 

Clearly,  each  internal  red-green  intersection  is  charged  to  at  most  one  cell.  Moreover, 
it  is  easily  checked  that  each  internal  red-green  or  red-blue-green  cell  is  charged  exactly 
C^  red-green  intersections,  namely  those  between  the  red  lines  whose  right  endpoint  lies 
in  Q'j  and  the  green  lines  whose  right  or  left  endpoint  (as  the  case  might  be)  lies  in  Qf 
(by  construction,  all  those  intersections  points  do  lie  inside  Q).   Hence  there  axe  at  most 

— ^  red-green  and  red-blue-green  cells.  Similarly  we  can  prove  that  there  axe  at  most  — ~- 

c  .  c^ 

red-blue  cells.  Thus,  the  overall  number  N  of  cells  after  the  third  step  is  bounded  by 


K, 


Krb 


N      <      2(*  +  u  +  t;)  +  -:^  +  -^-6 


^2  ^2 


<     2 


rrir 


+ 


rria 


+ 


rrii, 


+  ^^{Kr,  +  Kr,)  -  6 


f' 


<     -(m  +  2)  +  —{Krg  +  Krb)  -  6 
2  K 


<     -^m  +  --2 


Lemma  4.8  Let  M  he  the  total  number  of  cells  produced  by  the  above  algorithm.   Then 

.r  ^   3  2K 

C  C^  . 

Proof:  It  is  easy  to  see  that  for  any  two  adjacent  boundary  cells,  one  of  them  has  at  most 
four  edges,  therefore  there  are  at  most  \S\  cells  with  five  or  more  edges  (see  figure  10). 
Moreover  there  is  at  most  one  cell  with  seven  or  eight  edges.  Since  each  5  or  6-cell  is  par- 
titioned into  two  convex  quadrilaterals  and  the  7  or  8-cell  into  three  convex  quadrilaterals, 
the  fourth  step  creates  at  most  \S\  +  1  convex  quadrilaterals,  which  in  conjunction  with 
(4.7)  impUes  that  after  the  fourth  step  the  number  of  cells  is  at  most 


2m       a:      „      , 

-T-  +  -T^-2  +  (t  +  u  +  v-Z)   < 


3m      K 


Finally,  the  fifth  step  adds  at  most  [^^1  new  cells  at  each  quadrilateral  containing 
Ka  >  QVK  intersection  points.   Since  there  can  be  at  most  -A-  such  quadrilaterals,  the 
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Figure  10:  No  two  adjacent  cells  have  more  than  four  edges 


total  number  of  new  cells  is  at  most 

^IcVF  +  V  -    c    -  c^' 

because  we  apply  this  step  only  when  K  >  6^^.  Thus 

, ,       3m       2K 


D 


Next,  we  bound  the  maximum  number  of  lines  passing  through  any  cell  and  the  total 
number  of  line-cell  crossings. 

Lemma  4.9  For  each  1  <  i  <  M,  the  number  m,  of  lines  passing  through  the  cell  Q, 
satisfies 

m,   <  max{2V2\/K  +  4C,  2y/K +  &(;}.  (4.8) 

Proof:  Note  that  a  cell  Qi  produced  by  our  algorithm  is  a  subset  of  some  cell  .F^  G  F  (the 
cells  obtained  by  overlapping  red,  blue  and  green  pseudo  edges).  Therefore  it  suffices  to 
bound  the  number  of  lines  passing  through  a  cell  in  F.  First,  let  us  consider  a  red-green 
cell  ^ij  =  Qf  n  Q^.  By  Lemma  4.5,  Qf  meets  at  most  2JKgg  +  2^  green  lines  and  Q^  meets 
at  most  2\/Krr  +  2^  red  lines,  and  since  no  blue  hne  passes  through  a  red-green  cell,  there 
are  at  most  2{\/KrT  +  jKgg  -\-  2C,)  lines  pass  through  .F,j.  Similarly,  we  can  show  that  a 
red-blue  cell  meets  at  most  2{\/Ktt  +  y/K^  +  2C)  lines. 

Finally,  let  Tu  be  a  red-blue-green  cell,  then  Tu  is  the  intersection  of  Qj,  Q\  and  Q\. 
Since  Q\  (resp.  Q\)  meets  at  most  2^  green  (resp.  blue)  lines,  Tu  has  at  most  2\/Ktt  +  6C 
Unes.  Now  the  lemma  follows  from  the  fact  that  max{-yKV7  +  \/i^66,  >/^V7  +  J^gg]  ^ 
^/2K. 

U 

For  X,  y  G  {r,  6,  g],  let  cSf  denote  the  nvunber  of  pseudo  edges  of  color  x  intersected  by 
the  line  U  of  color  y.  It  follows  from  Lemma  4.2  that 


m. 


E^r  <  ^ 


+  mj. 


1=1 
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In  the  next  two  lemmas  we  bound  the  "mixed"  sums  XI^'i  ^I^,  HTJi  ^l''^  Y1T=\  ^f  ^-nd 

:=1  "i    ■ 


Lemma  4.10 


Y.^7  <^  +  ^9     and      Y.^1'  <   -^  +  rn,. 


t=l 


rfc+4 


(a)  (b) 

Figure  11:  Illustration  for  Lemma  4.10  and  Lemma  4.11 

Proof:  We  will  only  prove  the  first  paxt  of  the  lemma;  the  second  part  can  be  proved  in 
a  symmetric  way.  Let  ii  be  a  green  line  intersecting  (5[^  red  pseudo  edges,  r^+i,  . . .  ,  r^+^rsi, 
and  let  i'  be  a  red  line  having  its  right  endpoint  in  Q[,  for  k  +  \  <  I  <  k  +  8\^  (see 
figure  11a).  It  is  easily  checked  that  i'  intersects  £,  inside  Q.  Since  each  Q\  conteiins  the 
right  endpoints  of  C,  lines,  ii  intersects  at  least  {8^  —  1)C  hnes  inside  Q.  But  there  are  at 
most  Krg  red-green  intersection  points,  therefore  summing  over  all  green  lines,  we  obtain 

mg  Trig  T^ 

B-^r^-lX   <   Kr,     or      ^<5r    <    ^  +  m. 


«=1 


t=l 


c 


Lemma  4.11 


K, 


mr 


rg 


and      Yl^r'   ^ 


i=l 


Proof:  Let  £i  be  a  red  line  intersecting  (^f^  green  pseudo  edges,  gi,  ...  ,  gsar.  Suppose  that 
the  right  endpoint  of  £j  Hes  on  ^TgJ.  Let  £'  be  a  green  hne  whose  left  endpoint  (on  qiQ^)  lies 
in  Qf ,  for  /  <  ^f  (see  figure  lib).  It  is  easily  seen  that  i'  intersects  £,  inside  Q.  Similarly, 
if  the  right  endpoint  of  £,•  is  on  9394,  then  £,  intersects  any  green  line  whose  right  endpoint 
(on  9394)  lies  in  Qf,  for  /  <.■£?'".  Since  each  2f  contains  (  left  endpoints,  eind  also  (  right 
endpoints,  of  green  lines,  €,  intersects  at  least  6^  ■  (^  green  Unes  inside  Q.  Summing  over 
ail  red  lines,  we  obtain 


E^r-c  <  Kr,  or  Y.^r  < 


K, 


^9 


i=\ 


»=1 


c 


Similarly,  we  can  prove  that  EI^'i  ^Y  <  ^ 
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Lemma  4.12   Let  m\  denote  the  number  of  lines  passing  through  the  i      cell  ofF,  then 


N 


Y.<   <"^K+- 


m 


(4.9) 


i=l 


Proof:  Let  r/J  denote  the  number  of  cells  crossed  by  a  line  £_,  €  £x,  where  x  G  {r,  6,  g}. 
Each  cell  T  meeting  £_,  is  charged  to  the  leftmost  point  wjr  of  dJ^r\^j,  therefore  for  all  cells 
crossed  by  £_,,  except  the  leftmost  one,  wjr  is  an  intersection  point  of  a  pseudo  edge  and  ij 
(see  figure  12).  For  red  Unes  it  follows  from  Lemma  4.2  ajid  Lemma  4.11  that 


.'       C2 


>  ^ 

C4 


we-v^ 


W2    -  -'^       ^4     %/ 


C6 


Figure  12:  Line  ij  and  the  cells  in  F;  cell  c,  is  charged  to  w, 


E'?:  =  E(<5r  +  <5r  +  ^r  + 1) 


t=i 


i=l 


<     2^  +  ^  +  ^+2m, 
Similarly,  using  Lemma  4.2  and  Lenmia  4.10  we  can  show  that 


TTlft 

E 

i=l 


mft 


E^.*   =   E(^'*  +  <5[' +  ^f  + 1) 


«=i 


and      Er^f      <     2^  +  ^  +  3m, 
»=l  ^  s> 


N 

I 

:=1 


m. 


Since  ^m-  =  -^   E^J'  sunaming  (4.10),  (4.11)  and  (4.12),  we  obtain 


x&       j=l 
(r.b.g} 


N 


E  "^^     ^     7i^-r  +  Kbb  +  Kgg  +  Krg  +  i^.i)  +  2{mg  +  rrib)  +  2m, 


t=i 


<     fif+|m 


(  because  m,  >  y;  see  Section  2) 


(4.10) 


(4.11) 
(4.12) 
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Next,  we  bound  T.^i  "^^i  that  is  the  total  number  of  line-cell  crossings,  for  the  final 
cells  produced  by  the  algorithm. 

Lemma  4.13  Let  m,  denote  the  number  of  lines  in  cell  Qi;  then 

M  rr 

Y^  m.  <  (4  +  4v^)—  +  (5  +  4v^)m. 

Proof:  Let  rria  denote  the  number  of  Hues  meeting  the  a'''  cell  after  the  fourth  step.  First 
consider  the  case  when  F  has  a  cell  ^^  with  seven  or  eight  edges,  say  eight  edges  (see 
figure  8).  Let  ri_i  and  r,  denote  the  red  pseudo  lines  bounding  JT^.  In  this  case  the  fourth 
step  adds  only  two  pseudo  edges,  7£,  and  7h-  Since  jl  does  not  intersect  green  lines  and 
intersects  only  those  red  lines  whose  left  endpoints  lie  in  Q[,  it  adds  at  most  (  +  mi)  line-cell 
crossings.  Similarly  7/?  adds  at  most  (,'  -|-  rrig  line-cell  crossings.  Moreover,  if  F  has  a  7-cell, 
then  it  can  also  have  one  5-cell.  Thus  by  Lemma  4.12, 

v-<  2  5 

a  C  ^ 

Note  that  in  this  case  the  first  three  steps  create  at  least  three  pseudo  edges,  which  implies 
that  m  >  3C-  Since  nig  +  rrib  <  y,  we  have 

2 
Y'TTia   <   -K  +  5m. 

a  's 

On  the  other  hand  if  all  cells  after  the  third  step  have  at  most  six  edges,  then  every  cell 
is  partitioned  into  at  most  two  cells.  Thus  by  Lemma  4.12,  the  total  number  of  line-cell 
crossing  after  the  fourth  step  is 

4 
^ma  <  -K  +  5m. 

a  's 

The  final  step  creates  at  most  ^^  new  cells.    Moreover,  in  that  step  a  cell  J-a  is 

partitioned  only  when  K  >  6C^,  in  which  case  ma  <  2>/2\/^  +  4<^  (cf.  Lemma  4.9).  Hence, 
we  obtain 

M  A  2\/K 

^TTii     <     -K  +  5m  +  -^(2V2Vk  +  4C) 


t=i 


=     (4-h4v/2)y  +  sVk  +  5m 
<     (4  +  4^2)^  +  {5  +  4\/2)m 


because  K  <  ^ 
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Finally,  we  ajiaiyze  the  running  time  of  the  algorithm. 

Lemma  4.14  It  requires  O  Um  +  7-)log77ij  time  to  compute  all  subquadnlaterals  Qi,  . .  .  ,  Qi^, 
to  determine  the  lines  passing  through  each  Qi,  and  to  count  the  num,ber  of  intersection 
points  m  each  Qi. 

Proof:  It  follows  from  Lemma  4.5  that  Q"",  Q^  and  Q*"  can  be  computed  in  time  0{{m.r  + 
mb  +  m^)  log  m)  =  0{m  log  m).  Once  we  have  computed  all  pseudo  edges,  the  set  F  of  cells, 
obtained  by  overlaying  these  pseudo  edges,  can  be  computed,  for  example,  by  sweeping  a 
vertical  hne  through  the  quadrilateral  Q.  The  time  spent  in  the  sweep  is  0((|<5|  +  |F|)  log  5) 

=  O  ((—  +  — )logm  J,  because  it  encounters  at  most  0(|F|)  intersection  points.    Since 

step  4  simply  scans  all  boundaxy  cells  and  sphts  those  having  more  than  four  edges,  it  can 
be  easily  done  in  0(|5|)  time. 

As  for  the  final  step,  let  m^,  (resp.  Ka)  be  the  number  of  Hues  meeting  (resp.  intersection 
points  contained  in)  the  a*''  cell  produced  by  the  first  four  steps.  The  total  time  spent  in 
computing  the  hnes  passing  through  each  cell  J"^  is  at  most  0(m  log  m  +  XI^  m^  log  m)  = 
0{{m  +  -^)logm)  (of.  Lemma  4.13),  and  it  follows  from  Lemma  3.1  that  the  quantities 
Ka,  for  aU  cells  Ta,  can  be  calculated  in  time  0{{m  +  Y-)logm).  Moreover,  step  5  of 
the  algorithm  partitions  the  a"*  cell  into  at  most  f^l  convex  quadrilaterals,  which  is 
done  applying  the  algorithm  of  Section  3  at  most  f^]  -  1  times,  each  apphcation  taking 
0(mc  log  mc)  time  (cf.  Theorem  3.3).  Therefore,  the  total  time  spent  in  the  final  step  is 
bounded  by 

^(?(   CV^     -l)-^.logm.) 

(because  m„  <  2y/2\/K  +  4Q 
=    o(^ij  +  VK)\ogm) 

(K  \ 

(— +m)logmj 

(because  K  <  — ). 
Thus,  the  overall  running  time  of  our  algorithm  is  O  ((m  +  |^)logm). 
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Hence,  we  can  conclude  that 

Theorem  4.15  We  can  partition  the  convex  quadrilateral  Q  into  M  <  ^  +  ^  convex 
subquadrilaterals  Qi,  ...  ,  Qm,  with  the  property  that  each  Q,  is  crossed  by  m,  lines  and 
contains  Ki  intersection  points,  so  that  the  following  conditions  are  satisfied: 

(0     E^im,  <  Aim  +  A2f;     maxm.  <  Ag^F  +  A4C 
{li)     Ei^i  K,  =  K;  max  A'.  <  CV^ 


where  Ai  =  (5  +  A\/2),  Aj  =  (4  +  4^2),  A3  =  2v^  and  A4  =  6. 


D 


5      Partitioning  the  plane  into  quadrilaterals 

In  this  section  we  obtain  the  main  resiilt  of  this  paper.  Let  C  =  {^1,^2,  •  ■  ■  ,^n}  be  a  set 
of  n  lines  in  the  plane,  and  R  an  "enclosing  rectangle"  of  £,  i.e.  one  that  contains  aJl  (2) 
intersection  points  of  C.  For  a  given  integer  r  >  0,  we  want  to  partition  R  into  0{r^  log  r) 
convex  quadrilaterals,  for  some  constant  /?  >  0  to  be  determined  later,  so  that  each  of  them 
is  crossed  by  at  most  j  hnes. 

The  idea  is  to  use  the  algorithm  of  the  previous  section  recursively.  Fix  C  =  ^1  where 
7  is  a  constant  to  be  chosen  later  (for  simplicity  let  us  assume  that  n  is  a  multiple  of  jr). 
At  each  step  of  the  algorithm  we  have  a  convex  quadrilateral  Q,  which  meets  m  of  the 
given  Unes  and  contains  K  of  their  intersection  points.  Initially  Q  =  R,  m  =  n  and  K  = 
(2)  (assuming  that  the  lines  are  in  general  position).  At  each  recursive  step,  the  algorithm 
proceeds  as  follows:  If  m  <  ^,  there  is  nothing  to  do,  so  we  stop;  otherwise  partition  Q 
into  M  convex  quadrilaterals  using  the  algorithm  of  the  previous  section  (for  the  initial 
fixed  value  of  Q. 

Let  rrii  denote  the  number  of  lines  meeting  the  convex  quadrilateral  Q,-  and  let  Ki 
denote  the  number  of  intersection  points  of  £  contained  in  Q,.  Let  C(m,  K)  denote  the 
maximum  niimber  of  cells  into  which  such  a  Q  will  be  partitioned  by  all  subsequent  re- 
cursive applications  of  the  algorithm,  and  let  T{m,  K)  denote  the  maximum  time  required 
for  such  a  partitioning.  It  follows  from  Theorem  4.15  that 

M 


C{m,K)      <      < 


T{m,  k)     <      ' 


J2C{mi,Ki)     ifm>^ 
1  otherwise 

^T(mi,A:i)  +  D  fm  +  —  J  logm     if m  >  -  .^2) 


D  otherwise 
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where  D  is  some  constant  >  0  and  A/,  m,  and  A',  satisfy  the  bounds  given  in  Theorem  4.15. 

Next  we  bound  the  values  of  C{mJ\)  and  T{m,K)  using  (5.1)  and  (5.2).  In  what 
follows  by  Logx  we  mean  maxjloga:,  1},  and  similarly,  later  on,  by  LogLogx  we  mean 
max{loglogx,loglog  >/6}. 

Lemma  5.1    There  exists  a  constant  E  >  0  such  that 


C{m,K)  <E\-^  +  '-^\  Log'^^,  (5.3) 

where  (3  =  max{log  Ai,  Iog(l  +  A2)]  <  3.33. 

Proof:  We  prove  the  above  inequality  by  induction  on  K.   If  we  choose  7  >  4(a/3  4-1), 
then  for  K  <  6C^,  after  applying  the  algorithm  once,  we  have 

maxm.      <      max{2V2v/^4- 4C,  2^^"  +  6C} 
=     max{(4V3  +  4)C,(2v^  +  6)C} 

=     4(n/3  +  1)C      <      - 

r 

and  therefore,  the  algorithm  stops.   By  Theorem  4.15  this  step  partitions  Q  into  at  most 
3(^4-^)  cells.  Thus,  if  we  choose  E  >  3,  then  for  K  <  6C^  (5.3)  holds  trivially. 

For  K  >  6C^  suppose  inductively  that  (5.3)  holds  for  all  K'  <  K.  Since  K,  <  C\fK 
and  a:  >  6C^  Ki  <  ^  ■  \fK  =  ^  <  K.  Therefore  by  induction  hypothesis,  (5.3)  holds 
for  all  Qi,  so  (5.1)  implies 

C{m,K)     <     f  ^fr^  +  ^^Log^^  (5.4) 

Since  Vz,  Ki  <  C,\/K,  we  have 


c     cv      c 


Log—  <  Log-—  =  maxjlog  J  — ,  1}. 
K  >  6C^  implies  that  log  yj^  >  1,  therefore  (5.4)  can  be  written  as 

M  1      M 


(By  Theorem  4.15) 
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Let  0  =  max{log.4i,  log(l  +  .42)},  then 
Therefore 


Remark  5.2:  The  actual  value  oi  0  depends  on  the  constants  appearing  in  Theorem  4.15(i). 
We  believe  that  these  constants  are  not  optimal,  so  ^  is  Ukely  to  be  smaller  than  log2(5  + 
4v/2)  =  3.33. 

Next,  we  bound  the  running  time  T{m^  K)  of  the  algorithm. 
Lemma  5.3   There  exists  a  constant  F  >  Q  such  that 

(k\  k  k 

m  +  — jlogm-Log^— -LogLog— .  (5.5) 

Proof:  Again  we  prove  the  inequality  by  induction  on  K.  In  Lemma  5.1  we  showed  that, 
for  K  <  6C^,  the  edgorithm  stops  after  one  step  of  recursion.  By  Theorem  4.15,  this  step 
requires  0{{m  +  7-)  log  m)  time,  therefore  if  we  choose  F  sufficiently  large,  then  (5.5)  holds 
trivially. 

For  K  >  6C,  suppose  inductively  that  (5.5)  holds  for  ail  K'  <  K.   Since  Ki  <  K  (cf. 
Lemma  5.1),  by  induction  hypothesis,  (5.5)  holds  for  all  Q,-,  so  (5.2)  implies 

Af  /  IV-    \  TV-  IV- 

T{m,K)     <     5^F(m.  +  ^)logm,.Log^^-LogLog^  + 


t=i 


^y     o     .         o    ^2         o     0^2 


Z?fm  +  — jlogm  (5.6) 

Note  that,  for  all  i,  Ki  <  (^vK,  therefore 

K  K  Ik  fW 

Log^--^LogLog--j  <  max{log^  W  — ,  1}  •  max{loglog  W  — ,  log  log  V^}. 
But  K  >  6C^,  therefore  it  is  easily  seen  that 

/  ly-  j  LT  j  TV-  /  TV- 

^^S'^V  72^°sLogW—  =  log^  J—  loglog  , /— .  (5.7) 

V  C  V  C  V  C  V  c 
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Substituting  it  in  (5.6),  we  obtain 

T{m,K)     <     f^rfm. +  yjlogm-log'^  J— -loglog  J^  +  £>(m  +  yjlogm 

<     (  E  -^.  +  ^  g  ^«  J  ^  ■  I^S"  ^  (l°S  log  -  -  1  j  log  m  + 


i=i 


D  (  m  +  —  )  log  m 


<     L4im  +  (A2  +  1)  — j  2;^ -log^— flog  log— -ij  log  m  + 


D  [Tn^ )  logm         (by  Theorem  4.15) 


K 

m  +  — 


Flog^  ^  •  loglog^  +  (^D  -Flo^  p 


logTTZ 


because  f3  =  max{log  Ai,log(A2  +  1)}.     Since  K   >  6C  we  have  Flog^  |^   >  2'^,  and 
therefore,  by  choosing  F  =  ^,  we  obtain 

T{m,K)     <     F(m  +  — jlogm-log^  — -loglog  — 


<     F  (m  + yj  logm-Log^^    LogLog^. 

D 

These  lemmas  imply  that 

Theorem  5.4   Given  a  set  C  of  n  lines  in  the  plane  and  a  parameter  I   <  r  <  n,   we 
can  partition   an   enclosing  rectangle  R  into  0{r^  \o^  r)   convex  quadrilaterals   in   time 

0(nr  log  nlog^r  log  log  r),  so  that  each  quadrilateral  meets  at  most  —  lines  of  C. 

r 

Proof:  Choose  C  =  ^^  where  7  =  4(v^  +  1).  At  the  top  level  of  the  recursion  m  =  n  and 


K  =  [1)  <^.  Therefore 


K  nV2      _   7^2 


^2  „2/-y2r2  2 

-T  =  — ; —  =  "rr.  Substituting  these  values  in  (5.3)  we  obtain 


m  n 

and  --  = 


C(n,nV2)     <     F[7r  +  ^)Log^^ 


=     0(r2log^r) 
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K  n^/2         -yr 

and  by  substituting  —  =  -r—. —  =  -—n  in  (5.5)  we  get 

(,  2n/7r         4 


2^2  ^2^2 


n  +  — n  j  log  n  •  Log'^— ^LogLog 
=     O  fnr  log  nlog'^  r  log  log  rj 

n 

Remark  5.5:  To  m^ike  the  notations  easier  to  follow,  we  shall  henceforth  replace  the  term 
log'^  r  ■  log  log  r  by  log*^  r,  for  some  fixed  constant  u;  slightly  larger  than  0.  Since  we  do  not 
know  what  the  best  value  for  /?  is,  this  convention  involves  no  real  loss  of  information. 

6      Constructing  approximate  levels 

In  this  section  we  describe  the  second  phase  of  our  algorithm  that  reduces  the  number 
of  triangles  from  0(r^  log'^r)  to  0{r^),  while  maintaining  the  property  that  each  triangle 
meets  at  most  O(^)  lines  of  C.  As  mentioned  in  the  introduction,  this  second  phase  is  not 
required  in  most  of  the  apphcations.  As  an  intermediate  step,  the  algorithm  constructs  an 
^-approximate  leveling  of  A{C)  (as  defined  in  Section  2),  with  0{r^  log'^r)  edges  in  total, 
from  the  partition  obtained  in  the  first  phase  of  the  algorithm.  Once  an  ^-approximate 
leveling  has  been  constructed,  we  proceed  in  the  same  way  as  Matousek  [Maa].  We  first 
describe  how  to  obteun  an  ^-approximate  leveling. 

Let  Q  denote  the  planax  map  induced  by  the  preceding  partition  of  the  enclosing 
rectangle  R.  We  assiime  that  all  lines  in  C  intersect  5R  at  its  vertical  edges.  Let  A 
=  { 01,02,  ...  ,an}  (resp.  B  =  {61,62?  •••  ^K})  denote  the  intersection  points  of  the  lines 
in  C  and  the  left  (resp.  right)  vertical  edge  of  5R,  sorted  in  decreasing  order  of  their  y- 
coordinates.  For  1  <  i  <  r,  add  Oin/r  and  6,„/r  to  the  set  of  vertices  of  Q.  We  triangulate 
all  faces  of  Q,  in  time  0{r^  log^  r).  The  triangvilated  map  Q  also  has  0{r'^  log^  r)  edges. 
The  following  observations  enable  us  to  compute  approximate  levels  efiiciently. 


Figure  13:  A  triangle  intersecting  both  ICa  and  ICi,- 
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Lemma  6.1  Lei  )Ca  (resp.  K,b)  denote  the  k  -  ^  (resp.  k  +  ^)  level  of  A(C).  Then  the 
interior  of  a  triangle  in  Q  does  not  intersect  both  Ka  arid  Kb- 

Proof:  Suppose  to  the  contrary  that  there  is  such  a  triangle  A  =  ViV^Vk  that  crosses  Ka 
and  Kb.  Let  p^  (resp.  pb)  denote  a  point  in  ^^  H  A  (resp.  Kb  fl  A)  (see  figure  13).  The 
segment  p^  hes  entirely  in  A,  but  by  Lemma  2.1  it  intersects  at  least  7  +  1  lines  of  £, 
which  contradicts  the  property  that  no  triangle  in  Q  intersects  more  than  7  lines  of  C. 


Lemma  6.2  Let  Ka  (resp.  Kb)  denote  the  k  -  ^  (resp.  A:  4-  ^)  level  of  A{C).  Let  vq  be  a 
vertex  of  Q  on  the  left  vertical  edge  of  K,  whose  level  is  between  k  -  f^  and  k  +  f^.  There 
exists  a  path  H  in  Q  from  Vq  to  the  right  vertical  edge  of  R  which  lies  between  Ka  and  Kb- 
Moreover  IT  can  be  converted  into  a  monotone  path  without  increasing  its  number  of  edges, 
such  that  the  new  path  also  lies  between  Ka  o,nd  Kb- 

Proof:  Let  A  =  Ai,  Aj ,  ...  ,  At  denote  the  sequence  of  triangles  visited  if  we  follow  Kb 
from  left  to  right  (if  a  portion  of  Kb  hes  on  an  edge  of  a  triajigle,  then  we  pick  up  the 
triangle  lying  below  that  edge).  If  a  triangle  appears  more  than  once  in  A,  then  retain 
only  its  first  occurrence.  It  is  easily  seen  that  A  forms  a  connected  region  from  the  left 
vertical  edge  to  the  right  vertical  edge  of  5R,  and  its  boundary  is  formed  by  the  edges  of  Q 
(see  figure  14);  we  call  it  a  corridor.  Let  11  denote  the  top  portion  of  5A.  11  is  a  connected 
polygonal  path  from  the  left  vertical  edge  of  5R  to  its  right  vertical  edge,  formed  by  the 
edges  of  Q.  Obviously  11  does  not  intersect  Kb-  Since  each  triangle  of  A  intersects  Kb,  11 
cannot  intersect  Ka  (cf.  Lemma  6.1),  thus  11  hes  between  Ka  and  Kb- 


Figure  14:  Level  Kb  contained  in  the  corridor  A;  IT  is  the  upper  boundary  of  A. 


To  prove  the  second  part,  let  T  denote  any  polygonal  path  lying  between  Ka  and  Kb- 
Let  vi  be  the  first  vertex  on  F  at  which  T  turns  backwards.  Let  ^  denote  the  first  point  on 
r  after  vi  (along  F),  which  has  the  same  i-coordinate  as  ui.  Since  v-i  lies  above  Kb,  and 
Kb  is  a  i-monotone  path,  the  vertical  segment  vi^  does  not  cross  Kb-  The  same  argument 
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implies  that  Vi^  does  not  intersect  /Ca,  and  therefore  the  new  path  also  hes  between  K,a 
and  fCh-  Whenever  we  add  a  new  edge,  we  remove  at  least  one  edge  of  F,  which  proves 
that  shortcuting  F  does  not  increase  its  number  of  edges. 

D 

Since  the  path  11  hes  between  the  k—^  and  k  +  ^  levels  o{A{C),  it  is  an  ^-approximate 
fc-level  of  A{C).  Therefore  Q  contains  an  ^-approximate  ^-level,  for  each  i  <  r.  Let  Qk 
denote  the  subgraph  of  Q  consisting  of  all  edges  that  fully  lie  between  the  levels  k  ~  j^  and 
k  +  Y-  '^°  obtain  Qk,  we  need  to  compute  the  level  of  each  vertex  of  Q  ajid  the  highest  and 
the  lowest  levels  crossed  by  each  edge  of  Q.  Let  rl;{v)  denote  the  level  of  a  vertex  v  E  Q. 
Let  vivj  denote  an  edge  of  Q  with  x(u,)  <  x(u_,).  If  there  cire  h  lines  of  C  intersecting  vTU] 
and  k  of  these  Unes  lie  above  u,-,  then 

iP{vj)  =  ^(vi)  -  k  +  (h  -  k)  =  rp{v,)-2k  +  h  (6.1) 

For  each  vertex  of  Q  lying  on  the  left  verticeJ  edge  of  (9R,  we  can  compute  its  level  by 
counting  the  number  of  lines  of  C  lying  above  it.  The  partitioning  algorithm  produces  the 
subset  of  lines  of  C  that  crosses  each  edge  of  ^,  so  it  is  trivial  to  count  how  many  of  them 
lie  above  the  left  endpoint  of  e.  The  levels  of  all  vertices  of  Q  are  now  easy  to  determine  by 
propagating  levels  from  left  to  right  along  the  edges  of  Q,  using  (6.1).  As  for  determining 
the  levels  crossed  by  an  edge  e  of  Q,  we  sort  the  lines  intersecting  e  along  the  edge.  Once 
we  know  the  level  of  the  endpoints  of  e,  we  can  easily  compute  the  levels  crossed  by  e. 
Next  we  partition  G  into  Qi,  ...  ,Qr,  and  find  a  path  from  a.n/^  to  bin/r  in  Gi  using  a  depth 
first  search  from  a,„/r,  for  i  <  r.  Finally,  if  any  of  the  resulting  paths  is  not  x-monotone, 
we  malce  it  monotone  by  shortcuting  all  edges  that  turn  backwards.  Since  the  edges  in 
Qi,  for  i  <  r,  axe  pairwise  disjoint,  there  are  at  most  O(r^log'^r)  edges  in  the  resulting 
^-approximate  leveling  of  A{C). 

The  correctness  of  the  algorithm  follows  immediately  from  the  above  discussion,  so  we 
only  have  to  analyze  the  running  time  of  the  algorithm. 

Lemma  6.3  Given  a  set  C  of  n  lines,  we  can  construct,  in  time  0{nr  log  n  log'^  r),  an 
^■approximate  leveling  of  A{C)  having  only  O(r^log^r)  edges  in  total. 

Proof:  By  Theorem  5.4,  the  planar  graph  Q  can  be  constructed  in  time  0{nr  log  n  •  log'^  r) 
and  it  has  only  0{r'^  log^  r)  triangles.  It  follows  from  the  above  discussion  that  it  takes 
0{nr  logn  log"^  r)  time  to  compute  the  level  of  each  vertex  and  the  levels  crossed  by  each 
edge  of  Q.  Therefore,  we  can  obtain  Qi,  for  z  <  r,  in  0(nr  log  n  log'^  r)  time.  The  depth 
first  search  taJces  only  0(r^  log^r)  time,  and  it  takes  the  same  amount  of  time  to  convert 
the  computed  paths  into  monotone  paths.  Hence,  the  lemma  follows. 

D 
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Remark  6.4:  Matousek  [Maa]  also  constructs  an  approximate  leveling  as  an  intermediate 
step,  but  in  a  direct  and  much  simpler  (albeit  inefficient)  manner.  His  algorithm  works 
roughly  as  follows. 

Partition  R  into  r^  vertical  slabs  each  containing  at  most  0{  ^ )  intersection  points  of 
£.  Let  Vo,  . .  •  ,  K-2  denote  the  vertical  edges  of  these  slabs,  and  let  a,,i,  . . .  ,  a,„  denote  the 
intersection  points  of  VJ  and  C  sorted  in  decreasing  order  of  their  y-coordinates.  For  each 
j  <  ^  and  0  <  i  <  r^  connect  ai,2jn/r  to  a,+i,2jn/r-  These  polygonal  paths  are  shown  to 
form  an  ^-approximate  leveling  of  A{C). 

The  problem  with  this  approach  is  that  the  approximate  leveling  heis  0{r^)  edges  in 
total,  and  the  time  needed  to  obtain  it  is  0{nr'^  log^  r),  which  is  substantially  dominated 
by  the  peirtitioning  of  R  into  vertical  slabs.  Partitioning  R  into  r^  slabs  is  done  to  ensure 
that  no  segment  ai^2jn/rO-i+i,2jn/r  crosses  too  many  levels.  Using  our  improved  partitioning 
technique  we  are  able  to  obtain  an  ^-approximate  partitioning  that  has  almost  an  order 
of  magnitude  fewer  edges  (in  terms  of  r),  and  the  running  time  of  our  algorithm  is  also 
about  an  order  of  magnitude  faster.  If  r  is  small,  e.g.  0(1),  then  Matousek  algorithm  is 
better  (it  runs  in  optimal  linear  time),  but  for  large  values  of  r  it  becomes  very  inefficient. 
As  we  will  see  in  [Aga],  in  most  of  the  applications  it  is  desirable  to  use  a  large  value  of  r. 

After  computing  a  ^-approximate  leveling  of  A{C),  we  apply  the  same  technique  of 
Matousek  to  partition  R  into  0{r^)  triangles.  Let  Ki,K2,  ...  ,Kr  denote  the  set  of  ^- 
approximate  levels  of  A{C).  Using  Lemma  2.2  and  Lemma  2.3,  Matousek  proved  that 

Lemma  6.5  ([Maa])    There  exists  a  polygonal  path  H,  between  /f3,_i  and  K^i+i  such  that 


R" 


K.2 
K,3 


■  ■■  .       '      .-■•■■ 

•■     1  ■         1 

'■•••.    I 
1-J .J 

..    '1- 1 ^T 1 1 

--I .  ■  ■  ■  1  ■ -,1 

Figure  15:  Triangulation  of  /C* 

Note  that  K^i-x  Ues  between  the  levels  3i^  -  |  •  ^  and  Sz^  -|- 1  •  7,  and  /Ca.+i  Hes  between 
the  levels  3i^-i-^,  3i^  +  |-^.  Therefore  H.  Ues  between  the  kvels  3i^  q:  §  •  ^,  which 
shows  that  IIi,  ...  ,11^/3  forms  an  |^- approximate  leveUng  of  A{C).    By  applying  Suri's 
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jjgorithm  [Su]  of  computing  a  minimum  link  path  in  the  simply  connected  region,  lying 
between  A's,-!  and  A'si+i,  we  obtain  a  path  11',  such  that  in^l  <  |n,|.  Since  H',  lies  between 
iv3,_i  and  iv3i+i,  it  is  also  a  |^-approximate  ^ii^-level.  Hence,  we  have 

Corollary  6.6  Given  a  set  C  of  n  lines  and  a  ■parameter  \  <  r  <  n,  we  can  compute  a 
set  of  ^-approximate  levels  of  A{C),  fC  =  {Il[,  ...  ,n'^/3}  in  time  0(nr  log  n  log"' r)  with 

the  property  that  E^=^i  |n-|  =  Oir^). 

a 

Matousek  has  also  proved  that 
Lemma  6.7  ([Maa])    There  are  at  most  0{nr)  intersections  between  K,  and  C. 


In  view  of  Lemma  6.7,  K  can  be  decomposed  into  0{r'^)  edges,  each  intersecting  at 
most  O(^)  lines  of  £.  Next,  construct  the  vertical  decomposition  K.*  of  K  by  drawing  a 
vertical  line  from  every  vertex  of  K  in  both  directions  until  it  meets  an  edge  of  /C  or  R 
(see  figure  15).  Since  every  vertical  edge  added  to  AC  U  R  lies  within  —  levels,  it  intersects 
at  most  ^  lines  of  C  Therefore,  every  trapezoid  of  K*  intersects  at  most  O(-)  lines  of  £. 
Finally,  partition  each  trapezoidal  cell  of  K*  into  two  triangles.  Hence,  we  can  obtain  the 
main  and  final  result  of  the  paper. 

Theorem  6.8  Given  a  set  C  of  n  lines  in  the  plane,  and  a  parameter  \  <  r  <  n,  we  can 
decompose  the  enclosing  rectangle  R  into  0{r^)  triangles  in  time  0(nr log nlog'^r),  for 
some  constant  u  <  3.33,  so  that  no  triangle  meets  more  than  O(-)  lines  of  C. 

D 

7     Coping  with  Degenerate  Cases 

In  this  section  we  show  how  to  modify  our  partitioning  algorithm  so  that  it  also  works 
in  degenerate  cases,  when  more  than  two  lines  of  £  axe  concurrent,  or  more  than  one 
-  intersection  point  of  lines  in  £  he  on  the  same  vertical  hne.  To  facilitate  these  modifications, 
we  first  need  to  redefine  some  of  the  terminology,  introduced  in  Section  2  and  Section  4. 

(i)  If  <  >  2  lines  of  £  pass  through  a  point  p,  then  we  consider  p  as  (j)  intersection 
points.  Therefore,  although  the  number  of  distinct  intersection  points  in  £  can  be 
less  than  {^j ,  the  sum  of  their  "weights" ,  as  just  defined,  is  still  h)  ■ 
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(ii)  The  level  of  a  vertex  of  ACH)  can  no  longer  be  uniquely  defined,  so  we  redefine  the 
fc-level  o{  ACH)  to  be  the  polygonal  path  formed  by  the  closure  of  the  (open)  edges 
of  A{T-i)  whose  level  is  k.  Let  po,  ■■■  ,Pm  be  the  vertices  of  a  A;-level  of  AiH).  Let  w, 
denote  the  weight  of  Pi,  and  for  any  6  <  m,  let  p,-^  be  the  vertex  of  a  fc-level  such  that 
Yl\'-o^  Wt  <  jS  <  X]f=o^<-  We  now  define  the  S-simplified  k-level  to  be  the  polygonal 
path  connecting  po  to  p.-,,  pi,  to  p.-^,  ...  ,Pi,  to  p^,  for  s  =  L^^^J ,  concatenated 
with  the  left  and  right  rays  of  the  A;-level  incident  to  po  and  pm  respectively. 

(iii)  For  a  convex  quadrilateral  Q,  we  use  C  to  denote  the  set  of  lines  passing  through 
the  interior  of  Q,  and  we  let  K  denote  the  total  weight  of  the  intersection  points 
contained  in  the  interior  of  Q. 

(iv)  The  set  A  (resp.  B)  of  left  (resp.  right)  endpoints  of  the  hues  in  £,  defined  in  Section  2, 
now  becomes  a  multi-set,  because  many  lines  can  have  a  common  endpoint.  If  two 
lines  £i  and  ij  have  a  common  left  (resp.  right)  endpoint  x  and  ii  lies  counterclockwise 
(resp.  clockwise)  from  ^j,  when  directed  from  x  into  Q,  then  7r(z)  <  7r(j)  (resp. 
a{i)  >  cr(j)).  In  other  words,  A  and  B  axe  ordered  in  the  way  they  should  be,  if  we 
shrink  Q  slightly. 

(v)  Finally,  the  quantity  6i  defined  in  Section  4  denotes  the  number  of  pseudo  edges 
whose  relative  interior  intersects  £,.  Similarly  we  define  the  quantities  6^^ . 

Next  we  briefly  describe  the  modifications  in  ovu"  algorithm  and  in  its  analysis,  required 
to  make  them  work  in  degenerate  cases  as  well;  we  leave  it  for  the  reader  to  fill  in  the  details. 

Observe  that  in  Section  3.1  we  actually  covmt  the  number  of  inversions  to  compute 
Krr.Kbr  and  Kgr,  which  gives  the  total  weight  of  the  intersection  points  (not  the  number 
of  distinct  intersection  points)  contained  in  Q,  so  this  paxt  of  the  algorithm  does  not  require 
any  modification.  However,  now  it  is  not  always  possible  to  find  a  vertical  fine  that,  for  a 
given  k  <  K,  has  exactly  k  intersection  points  in  Q  to  its  left  (of.  Section  3.2).  Instead, 
we  find  the  rightmost  vertical  Une  having  <  k  intersections  in  Q  to  its  left.  If  we  interpret 
the  order  of  the  Hnes  of  C  along  a  vertical  line  A  as  the  order  that  would  result  by  moving 
A  slightly  to  the  left,  then  it  is  easily  checked  that  the  procedure  described  in  Section  3.2, 
with  obvious  and  trivial  modifications,  would  produce  the  desired  line.  We  now  apply 
this  procedure,  in  the  final  step  of  our  general  partitioning  algorithm  (cf.  Section  4.2),  to 
partition  a  cell  ^a  into  subcells,  each  conteiining  <  C,\/K  intersection  point  in  its  interior. 

The  algorithm  of  Section  4.1  remains  the  same.  Since  many  elements  in  A  or  5  can  have 
the  same  value,  two  or  more  pseudo  edges  can  have  a  common  endpoint  (see  figvire  16a), 
or  even  completely  overlap;  however,  we  will  regard  overlapping  pseudo  edges  as  a  single 
edge.  It  can  be  checked  that  our  new  conventions  regarding  A  and  B  imply  that  a  line 
ii  ^  C  intersects  at  least  |7r(z)  -  a{i)\  lines  inside  Q,  that  is.  Lemma  4.1  still  holds,  which 
in  conjunction  with  the  new  definitions  of  K  and  8  implies  that  Lemma  4.2  also  continues 
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Figure  16:  Cells  in  degenerate  cases:  (a)  special  case;  (b)  general  case 

to  hold.  Lemma  4.4  is  trickier  to  adjust.  We  redefine  pi,  as  the  number  of  lines  passing 
through  Qi  and  not  having  an  endpoint  whose  rank  in  A  or  in  S  is  between  (z  —  1)^  +  1 
and  iC-  With  some  caxe,  the  proof  can  be  modified  to  yield  the  same  bound  on  fii. 

Our  general  algorithm  described  in  Section  4.2  also  does  not  change  except  that  in  the 
final  step  we  now  use  the  modified  algorithm  mentioned  above  to  partition  the  cells  J^^ 
into  [w^l  cells,  each  containing  at  most  C,yK  intersection  points  in  its  interior.  Note 
that  now  we  can  have  some  triangular  ceils  that  do  not  have  an  edge  lying  along  dQ  but 
only  have  a  vertex  lying  on  dQ\  these  cells  are  also  considered  as  boundary  cells.  We  may 
even  have  some  zero-axea  cells,  caused  by  overlapping  pseudo  edges,  but  we  discard  these 
cells.  In  view  of  our  convention  regarding  the  weights  of  intersection  points  and  the  order 
of  the  lines  along  5Q,  it  is  easy  to  see  that  the  charging  scheme  of  Lemma  4.7  still  works. 
(More  specifically,  we  charge  a  red-green  intersection  to  a  red- green  cell  !Fi]  if  the  rank 
of  the  right  endpoint  of  the  red  line  is  between  (j  —  1)C  +  1  and  jC,  and  the  rank  of  the 
appropriate  endpoint  of  the  green  line  is  between  (i  —  1)(^-|-1  and  iC,,  and  similarly  for 
red-blue  intersections.)  Moreover,  using  the  same  argument  a^  in  Section  4.3  we  can  prove 
Lemma  4.10  and  Lemma  4.11  in  this  degenerate  setting  because  we  are  not  counting  those 
pseudo  lines  that  intersect  £  at  their  endpoints.  Finedly,  it  can  be  shown  that  Lemma  4.9, 
Lemma  4.12  and  Lemma  4.13  also  hold,  because,  although  a  line  £,  can  meet  the  boundary 
of  many  cells,  it  meets  the  interior  of  exactly  6i  -\-  1  cells.  Lemma  4.14  can  also  be  shown 
to  hold,  using  our  notational  convention  in  degenerate  cases. 

As  for  the  second  phase  of  our  algorithm,  Lemma  6.1  and  Lemma  6.2  are  not  aifected  by 
degeneracies  in  C.  However,  computing  the  graphs  Qk  now  becomes  slightly  more  difficult 
because  the  level  of  a  vertex  u  of  ^  may  be  undefined  (that  is,  when  u  is  a  "heavy- weighted" 
vertex  of  A{C)  as  well).  But  for  each  vertex  u  of  ^,  we  can  compute  the  lines  of  £  passing 
through  it  in  time  0{nr  log  n),  because  a  line  passing  through  a  vertex  v  either  lies  in  the 
interior  of  a  triangle  adjacent  to  v,  or  contains  on  of  the  edges  adjacent  to  v.  Moreover 
for  each  edge  incident  to  u,  we  can  compute,  in  time  0(nr  log  n),  its  level  in  a  sufficiently 
small  neighborhood  of  v.    Therefore,  we  can  still  propagate  the  levels  from  left  to  right 
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along  the  edges  of  Q,  and  can  determine  the  levels  crossed  by  each  edge  of  Q.  Thus,  we 
can  partition  Q  into  Qi,  ...  ,Qt  and  obtain  an  ^-approximate  levehng  of  A(C).  It  can  be 
checked  that  the  technical  results  in  [EW]  and  [Maa]  can  be  extended  to  the  degenerate 
case,  which  in  turn  implies  that  Lemma  2.3  and  Lemma  6.5  still  hold  if  we  follow  our  new 
conventions  regarding  simplified  levels  and  the  weights  of  intersection  points.  It  is  also 
easy  to  check  that  the  proof  of  Lemma  6.7  does  not  require  the  lines  of  £  to  be  in  general 
position. 

Hence,  we  can  conclude  that  Theorem  6.8  holds,  with  appropriate  modifications  as 
discussed  above,  even  if  the  lines  axe  not  in  general  position. 

8      Conclusions 

In  this  paper  we  presented  a  deterministic  algorithm  that,  given  a  set  £  of  n  lines  and  a 
parameter  1  <  r  <  n,  partitions  the  plane  into  0{r^)  triangles,  each  of  which  meets  at 
most  0{j)  lines  of  £.  Although  we  showed  that  our  algorithm  is  optimal  up  to  a  poly-log 
factor,  there  are  some  questions  that  stiU  remain  unanswered: 

(i)  The  time  complexity  of  our  algorithm  has  an  extra  log'^  r  factor.  The  value  of  u 
depends  on  the  constants  appearing  in  the  bounds  of  max  m^  and  J2tti  "^t  in  The- 
orem 4.15.  We  believe  the  values  of  these  constants  can  be  improved  by  doing  a 
more  careful  analysis,  so  a  natural  question  is  to  determine  the  best  possible  values 
of  these  constants. 

(ii)  As  mentioned  in  the  introduction,  we  conjecture  that  the  second  phase  of  our  al- 
gorithm is  redundant,  i.e.  the  first  phase  itself  or  some  appropriate  variant  of  it 
produces  0{r'^)  triangles;  so  far  we  have  not  been  able  to  prove  this  conjecture. 

(iii)  One  drawback  of  our  algorithm  is  that  its  space  complexity  is  0{nr).  This  is  because 
at  every  level  of  recursion  we  need  to  determine  the  fines  passing  through  each  cell. 
Can  the  space  complexity  be  reduced  to  0{n  +  r^)? 

(iv)  The  Q(nr)  lower  bound  applies  only  if  we  want  to  report  the  lines  passing  through 
each  triangle.  No  non-trivial  lower  bound  is  known  if  these  crossing  need  not  be 
reported.  A  challenging  open  problem  is  to  establish  a  similar  lower  bound  in  this 
case  too. 

Besides  the  above  open  problems  related  to  our  algorithm,  there  are  several  other  open 
problems  related  to  the  general  partitioning  problem. 

(i)  The  most  challenging  open  problem  is  to  generalize  our  algorithm  to  three  dimen- 
sions.   That  is,  for  a  given  set  H  oi  n  planes  in  R^  and  a  parameter  1  <  r  <  n. 
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how  fast  can  one  partition  space  into  0{r^)  tetrahedra  so  that  no  tetrahedron  meets 
more  than  O(-)  planes  of  H?  The  algorithm  of  [CF]  yields  such  a  partitioning  in 
time  0{n^°).  Since  the  number  of  plane-tetrahedron  crossings  is  Q{nr'^)  in  the  worst 
case,  an  idccd  solution  would  be  an  algorithm  whose  running  time  is  close  to  0{nr^). 

(ii)  Our  algorithm  relies  on  the  fact  that  the  objects  axe  straight  lines,  so  it  is  still  an 
open  problem  whether  a  similax  algorithm  exists  for  curves.  Note  that  the  algorithm 
of  Chazelle  and  Friedman  [CF]  works  even  for  curves,  but  their  construction  is  very 
inefficient.  Moreover,  in  the  case  of  curves,  the  techniques  of  [CF],  as  well  as  the  ran- 
dom sampling  technique,  both  give  a  slightly  weaker  bound,  namely,  they  partition 
the  plane  into  0{r^)  simply-shaped  regions,  each  meeting  at  most  0{j\ogr)  curves 
(instead  of  0(^)). 

(iii)  Suppose  we  have  a  collection  of  m  blue  lines  and  a  collection  of  n  red  lines.  For 
a  given  parameter  r,  how  feist  can  one  partition  the  plane  into  0{r'^)  triangles,  so 
that  each  tricingle  meets  O(^)  blue  lines  and  0(7)  red  lines?  Note  that  the  random 
sampling  technique  easily  yields  a  similax  partitioning,  with  each  triangle  meeting 
O(^logr)  blue  lines  and  0{jlogr)  red  lines. 

(iv)  Suppose  we  have  a  collection  Q  oi  n  segments  which  intersect  at  K  points.  The 
random  sampling  technique  shows  that,  for  a  given  parameter  r,  we  can  partition 
the  plane  into  0{r  -f-  ^^)  triangles  so  that  no  triangle  meets  more  than  O(-logr) 
segments  of  Q.  Can  our  algorithms  be  extended  to  yield  a  deterministic  algorithm 
that  produces  that  many  triangles,  each  meeting  0{j)  segments,  and  runs  in  time 
close  to  0{n  -\-  ^)  (the  lower  bound  on  triangle-segment  crossing  in  this  case  is 
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